ON CLASSIFICATION OF MODULAR TENSOR CATEGORIES 
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Abstract. We classify all unitary modular tensor categories (UMTCs) of rank 
< 4. There are a total of 35 UMTCs of rank < 4 up to ribbon tensor equivalence. 
Since the distinction between the modular S'-matrix S and —S has both topolog- 
ical and physical significance, so in our convention there are a total of 70 UMTCs 
of rank < 4. In particular, there are two trivial UMTCs with S = (±1). Each 
such UMTC can be obtained from 10 non-trivial prime UMTCs by direct prod- 
uct, and some symmetry operations. Explicit data of the 10 non-trivial prime 
UMTCs are given in Section 5. Relevance of UMTCs to topological quantum 
computation and various conjectures are given in Section 6. 



1. Introduction 

A modular tensor category (MTC) in the sense of V. Turaev determines uniquely 
a (2+l)-topological quantum field theory (TQFT) [Tuj (a seemingly different defi- 
nition appeared in [MSI] .) The classification of MTCs is motivated by the applica- 
tion of MTCs to topological quantum computing [F] [KH] [FKW] [FLWlj jFKLWj [P] , 
and by the use of MTCs in developing a physical theory of topological phases of 
matter [Wflj [MR] fFWTW] [Ki2] jWa] pVe] [DFNSS] . G. Moore and N. Seiberg ar- 
ticulated the viewpoint that rational conformal field theory (RCFT) should be 
treated as a generalization of group theory [MS2j . The algebraic content of both 
RCFTs and TQFTs is encoded by MTCs. Although two seemingly different def- 
initions of MTCs were used in the two contexts [MSlJ|Tu|, the two notions are 
essentially equivalent: an MTC in [MST] consists of essentially the basic data of 
a TQFT in [Wal] . The theory of MTCs encompasses the most salient feature of 
quantum mechanics in the tensor product: superposition. Therefore, even with- 
out any applications in mind, the classification of MTCs could be pursued as a 
quantum generalization of the classification of finite groups. 

There are two natural ways to organize MTCs: one by fixing a pair (67, A), where 
6? is a compact Lie group, and A a cohomology class G H 4 (BG; Z); and the other 
by fixing the rank of an MTC, i.e. the number of isomorphism classes of simple 
objects. If a conjecture of E. Witten were true, then every MTC would come from a 
Chern-Simons-Witten (CSW) TQFT labelled by a pair (G, A) [Witt] [M5T] [HKW] . 
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Classification by fixing a compact Lie group G has been carried out successfully for 
G=finite groups JIM] [FQ] , G = T n torus [Ma] [EM], and G = A,B,C,D simple 
Lie groups [FKJ [KWJ |TW] . In this paper, we will pursue the classification by fixing 
the rank. This approach is inspired by the study of topological phases of matter 
and topological quantum computing. Another reason is that we have evidence 
that there might be exotic (2+l)-TQFTs other than CSW theories jHRWj . 

Topological phases of matter are like artificial elements. The only known topo- 
logical phases of matter are fractional quantum Hall liquids: electron systems 
confined on a disk immersed in a strong perpendicular magnetic field at extremely 
low temperatures [Wilj [DFNSS] . Electrons in the disk, pictured classically as or- 
biting inside concentric annuli around the origin, organize themselves into some 
topological order [Wenj [WW lj |WW2] . Therefore, the classification of topological 
phases of matter resembles the periodic table of elements. The periodic table 
does not go on forever, and simpler elements are easier to find. The topological 
quantum computing project is to find MTCs in Nature, in particular those with 
non-abelian anyons. Therefore, it is important that we know the simplest MTCs 
in a certain sense because the chance for their existence is better. 

There is a hierarchy of structures on a tensor category: rigidity, pivotality, 
sphericity. We will always assume that our category is a fusion category: a rigid, 
semi-simple, C-linear monoidal category with finitely many isomorphism classes 
of simple objects, and the trivial object is simple. It has been conjectured that 
every fusion category has a pivotal structure [ENOj . Actually, it might be true 
that every fusion category is spherical. Another important structure on a tensor 
category is braiding. A tensor category with compatible pivotal and braiding 
structures is called ribbon. In our case a ribbon category is always pre-modular 
since we assume it is a fusion category. For each structure, we may study the 
classification problem. The classification of fusion categories by fixing the rank 
has been pursued in |01j|U2j . Since an MTC has considerably more structures 
than a fusion category, the classification is potentially easier, and we will see that 
this is indeed the case in Sections [3] and HI The advantage in the MTC classification 
is that we can work with the modular 5* matrix and T matrix to determine the 
possible fusion rules without first solving the pentagon and hexagon equations. 
For the classification of MTCs of a given rank, we could start with the infinitely 
many possible fusion rules, and then try to rule out most of the fusion rules by 
showing the pentagon equations have no solutions. However, pentagon equations 
are notoriously hard to solve, and we have no theories to practically determine 
when a solution exists for a particular set of fusion rules (Tarski's theorem on the 
decidability of the first-order theory of real numbers provides a logical solution). 
So being able to determine all possible fusion rules without solving the pentagon 
equations greatly simplifies the classification for MTCs. As shown in [HHj, all 
structures on an MTC can be formulated as polynomial equations over Z. Hence 
the classification of MTC is the same as counting points on certain algebraic 
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varieties up to equivalence. But all the data of an MTC can be presented over 
certain finite degree Galois extension of Q, probably over abelian Galois extension 
of Q if normalized appropriately. Therefore, the classification problem is closer 
to number theory than to algebraic geometry. The argument in Sections [3] and 
H] is basically Galois theory plus elementary yet complicated number theory. To 
complete the classification, we need to solve the pentagons and hexagons given 
the fusion rules. A significant complication comes from the choices of bases of 
the Horn spaces when solving the pentagon equations. The choices of basis make 
the normalization of 6j symbols into an art: so far no computer programs are 
available to solve pentagons with a fusion coefficient > 1, but one set of such 
fusion rules is solved completely [HHj . Currently, there are no theories to count 
the number of solutions of pentagon equations for a given set of fusion rules without 
solving the pentagons. For unitary MTCs, there is tension between two desirable 
normalizations for 6j symbols: to make the F matrix unitary, or to present all 
data of the theory in an abelian Galois extension of Q. For the Fibonacci theory, 
unitarity of the F matrix and abelianess of the Galois extension of Q cannot 
be achieved simultaneously, but with different F matrices, each can be obtained 
separately [FWj . This is the reason that we will only define the Galois group of a 
modular fusion rule and a modular data, but not the Galois group of an MTC. 

The main result of this paper is the classification of MTCs with rank=2, 3 and 
unitary MTCs of rank=4. The authors had obtained the classification of all unitary 
MTCs of rank < 4 in 2004 [Waj . The delay is related to the open finiteness 
conjecture: There are only finitely many equivalence classes of MTCs for any 
given rank. By Ocneanu rigidity the conjecture is equivalent to: There are only 
finitely many sets of fusion rules for MTCs of a given rank. Our classification 
of MTCs of rank < 4 supports the conjecture. We also listed all quantum group 
MTCs up to rank < 12 in Section Two well-known constructs of MTCs are the 
quantum group method, and the quantum double of spherical tensor categories 
or the Drinfeld center. The quantum double is natural for MTCs from subfactor 
theory using Ocneanu's asymptotic inclusions [EKj . It seems that this method 
might produce exotic MTCs in the sense of |HRWj . 

Our main technique is Galois theory. Galois theory was introduced into the 
study of RCFT by J. de Boer and J. Goeree |dBG| . who considered the Galois 
extension K of Q by adjoining all the eigenvalues of the fusion matrices. They 
made the deep observation that the Galois group of the extension K over Q is 
always abelian. This result was extended by A. Coste and T. Gannon who used 
their extension to study the classification of RCFTs [CGj . Fusion rules of an 
MTC are determined by the modular S"-matrix through Verlinde formulas. It 
follows that the Galois extension K is the same as adjoining to Q all entries of 
the modular S matrix. When a Galois group element applies to the S matrix 
entry-wise, this action is a multiplication of S by a signed permutation matrix, 
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which first appeared in |CG] . It follows that the entries of the S matrix are the 
same up to signs if they are in the same orbit of a Galois group element. For a 
given rank < 4, this allows us to determine all possible S'-matrices, therefore, all 
possible fusion rules. 

Note that the Galois group of a modular data does not change the fusion ma- 
trices, but it can change a unitary theory into a non-unitary theory. For example, 
the Galois conjugate of the Fibonacci theory is the Yang-Lee theory, which is non- 
unitary. We might expect that for each modular data, one of its Galois conjugates 
would be realized by a unitary MTC. This is actually false. For example, take a 

rank=2 modular data with S = [ \ \\, and T = [ „ . | . No Galois actions 

\-l -Ij \U ij 

can change the S matrix, hence the quantum dimension of the non-trivial simple 
object from —1 to 1, though the same fusion rules can be realized by a unitary 
theory: the semion theory. Reference [Rolj contains a set of fusion rules which 
has non-unitary MTC realizations, but has no unitary realizations at all. 

The paper is organized as follows. In Section [2J we study the implications of 
the Verlinde formulas using Galois theory. In Sections [3] and HJ we determine all 
self-dual modular S matrices of modular symbols of rank=2, 3, and unitary ones 
for rank=4. Rank=2 is known to experts, and rank=3 fusion rules have been 
previously classified |GPj . For modular data, Theorems 13.11 and 13.21 can also be 
deduced from [01] [02] . In Section we determine all UMTCs of rank < 4. In 
Section El we discuss some open questions about the structure and application of 
MTCs. In Appendix together with S. Belinschi, we determine all non-self dual 
unitary modular data of rank < 4. 

We summarize the classification of all rank < 4 unitary MTCs into Table [TJ 
There are a total of 70 unitary MTCs of rank < 4 (a total of 35 up to rib- 
bon tensor equivalence). The count is done in Section 15.41 Each such UMTC 
can be obtained from 10 non-trivial prime UMTCs by direct product, and some 
symmetry operations. The 10 non-trivial prime UMTCs are the semion MTC, 
the Fibonacci MTC or (Ai,3)i, the Z 3 MTC, the Ising MTC, the {A u 2) MTC, 
the even half of an SU(2) MTC at level 5 or (Ai,5)i, the Z 4 MTC, the toric 
code MTC, the (D 4 , 1) MTC, and the even half of an SU(2) MTC at level 7 
or (Al,7)i. Their explicit data are listed in Section IBTBl Out of the 10 non- 
trivial prime UMTCs, 9 are quantum group categories for a simple Lie group: the 
semion=S77(2)i, the Fibonacci=(G 2 )i, the Z 3 =577(3)i, the Ising=complex conju- 
gate of {E 8 ) 2 , the (A 1 ,2)=SU{2) 2 , the Z 4 = 5C/(4)i, the toric code= Spin{16) u 
the (-D 4 , 1) = Spin(8)x, and the (Ai, 7)i=complex conjugate of (^2)2- The Ising 
MTC and the SU(2) 2 MTC have the same fusion rules, but the Frobenius-Schur in- 
dicators of the non-abelian anyon a are +1,-1, respectively. The toric code MTC 
and the Spin(8)x MTC have the same fusion rules, but the twists are {1, 1, 1, — 1}, 
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Table 1. Unitary MTCs of rank< 4 





A 2 
1 










A 4 

Z 2 




N 4 
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Z 3 
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N 4 
(^i,5)i 

U 




A 10 

Z 2 x Z 2 


A 8 
Z 4 


N 8 
(At,3) 
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N 4 
(^i,7) i 

U 


N 6 
Fib x Fi6 

U 



and {1,-1,-1,-1}, respectively. We choose q = e « in the quantum group con- 
struction. In the Ising case, it is the q = e -2 ? theory for Eg at level=2. For 
notation and more details, see Section 15.31 We do not know how to construct 
(Ai, 5) i by cosets of quantum group categories. 

The information for each rank is contained in one row of Table [H Each box 
contains information of the MTCs with the same fusion rules. The center entry in 
a box denotes the realization of the fusion by a quantum group category or their 
products. We also use Fib to denote the Fibonacci category (Al,3)i. The upper 
left corner has either A or N, where A means that all anyons are abelian, and N 
that at least one type of anyons is non-abelian. The right upper corner has a num- 
ber which is the number of different unitary theories with that fusion rule. If the 
lower right corner has a U, it means that at least one type of anyons has universal 
braiding statistics for topological quantum computation. The detailed information 
about which anyon is abelian or non-abelian, universal or non-universal is given 
in Section 16.31 It is worth noticing that the list of all fusion rules up to rank=4 
agrees with the computer search for RCFTs in [GK] . We believe this continues to 
be true for rank=5. The rank= 6 list in [GK] is not complete. 

Finally, we comment on the physical realization of UMTCs. The existence of 
abelian anyons in v = ~ FQH liquids is established theoretically with experimental 
support, while non-abelian anyons are believed to exist at the v = | and v — y 
plateaus (see [DFNSS] and the references therein). Current experimental effort is 
focused on FQH liquids at v = |. But the fermionic nature of electrons complicates 
direct application of MTCs to FQH liquids because only anyonic properties of 
bosonic systems can be described fully by MTCs. In other words, we need a 
refined theory, e.g. a spin MTC, to describe a fermionic system |DWJ[BM]. 
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2. Galois theory of fusion rules 

In this section, we study the implication of Verlinde formulas for fusion rules of 
MTCs. For more related discussion, see the beautiful survey [G] . 

Definition 2.1. (1) A rank=n label set is a finite set £ of n elements with a 
distinguished element, denoted by 0, and an involution": L — * L such that 
= 0. A label i G £ is self dual if % = i. 

The charge conjugation matrix is the n x n matrix C = (<^). Note that 
C is symmetric and C 2 = I n , the n x n identity matrix. 

(2) A rank=n modular fusion rule is a pair (N;S), where N is a set of n 
n x n matrices iVj = (nf j)o<j,k<n-u indexed by a rank=n label set £, with 
n^j G Q , and S = (%)o<ij<n-i is an n x n matrix satisfying the following: 

(a) §oo = 1, s i = 5ij, and all s i)0 's are non-zero; 

(b) If we let D = JY2=o *?,o> then S = ^ is a symmetric, unitary matrix. 
Furthermore, the matrices iVj in N and S are related by the following: 

(2.1) N.S^SAi 

for all i G £, where A, = (S ab X ia ) nxn is diagonal, and \ ia = ff. 

The identities ( 12.1 p or equivalently the Verlinde formulas (I2.3P below 
imply many symmetries among •: rin ■ = 5jk, n^A = n^j = n~~. = n 3 . ~ . 

The matrix iVj will be called the ith fusion matrix. From identities (12. ip . 
the diagonal entries in Aj are the eigenvalues of Aj, and the columns of S 
are the corresponding eigenvectors. The non-zero number D will be called 
the total quantum order, di = Sio the quantum dimension of the i-th label, 
and D 2 = Yh=q d 2 the global quantum dimension. 

(3) A rank=n modular symbol consists of a triple (N; S,T). The pair (N; 

is a rank=n modular fusion rule with all n^- G N = {0, 1, 2, • • • } (here Soo 
is the (0,0)-entry of the unitary matrix S = (sr,-)o<ij<n-i), an d the n x n 
matrix T = (S a b9 a ) nxn is diagonal, and 9 = 1. Furthermore, S and T 
satisfy 

(i) (ST) 3 = (D + s 00 )S 2 ; 

(ii) S 2 = C; 

(iii) 6i G U(l) and = 6^ for each i; 

where D± = Yh=o ^i^ 2 ■ The following identity can be deduced: 

(2.2) D + D_ = D 2 . 

The complex number di will be called the twist of the iih. label. Note that 
Soo might he —j). A modular symbol is called unitary if each quantum 
dimension di is the Frobenius-Perron eigenvalue of the corresponding fusion 
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matrix iVj. In particular, the quantum dimensions dj's are positive real 
numbers > 1. 

(4) A modular symbol (N; 5*, T) is called a modular data if there is an MTC 
whose fusion rules, modular S"-matrix, and T-matrix are given by IN, T 
of the modular symbol. 

(5) Let A = {^ij}i,jeL for a rank=n modular fusion rule, and let K = Q(\ij), i, j G 
£j be the Galois extension of Q. Then the Galois group G of the Galois 
field K over Q is called the Galois group of the modular fusion rule. 

We are interested in searching for n+2 tuples (N Q , ■ ■ • , iV n _i; S, T) related in the 
correct fashion. We will index the rows and columns of matrices by 0, 1, • • • ,n — 1. 
Since Ni§ = SAi, the columns of S must be eigenvectors of Ni with eigenvalues X ij0 , 
X it x, • • • , and A i>n _i, respectively. Looking at the first entries of these columns and 
of NiS, and using the only non-zero 1 of the first row of Ni, we see that A^o = di, 
and djXij = §ij. It follows that K is the same as Q(%),i,j G £. Since S is 
symmetric, we see that for i ^ j we have djXij = diXj^, and sy = diXj^ = djXij 
for all i and j. Let n*„- denote the (j, k) entry of N{. Since 



(X i>0 o ••• \ 
Ai, a ••• 

V • • • A iin _J 



we compute for < j, k < n — 1 



f o\ „k _ \ " a i,m^,m^,m n -2 _ \ \ 1 m 

m=0 m=0 ^ 



The fusion matrices can also be described equivalently by fusion algebras. For a 
rank=n fusion rule, each label % is associated with a variable X{. Then the fusion 
ring R is the free abelian ring Z[X , ■ ■ ■ , X n _i] generated by X^s modulo relations 
(called fusion rules) XiXj = Yh=o n i,j-^ k - The fusion algebra will be F = R®z K, 
where K is the Galois field of the fusion rules above. We may replace K by C. If 
the modular fusion rule is realized by an MTC, then Xi is an equivalence class of 
simple objects, and the multiplication XiXj is just the tensor product. 

There are modular symbols that are not modular data. 

Example 2.2. Take the following 

s = i 

2 

and T = Diag(l,#, —1). The fusion matrices iVj are determined by the formulas 
( 12. 31) . hence are independent of 9. They are the same as those of the Ising MTC 





V2 









-y/2 




-V2 
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in Section IBTBl Therefore, for any 9 G U(l), we get a modular symbol. But only 
when 6 is a 16th root of unity, do we have modular data. 

Very likely the modular symbol of an MTC determines the MTC , and we do 
not know when a modular symbol becomes a modular data. 

Proposition 2.3. If(N;S, T) is a modular data, then we have: 

(1) OiOjSij = J2k n ij s ko9k; 

(2) {[,()■ ' oj^'^: 

where As* = 2n , ~n), + nLnV>; 

.02 

(3) Let Vk = Tyz Si,je£ n k jdidj-fe, then is if k ^ k, and is ±1 if k — k. Vk 

j 

is called the Frobenius-Schur indicator of k. 

(4) -D+Soo — e 2 ^ for some c G Q. The rational number c mod 8 is called the 
topological central charge of the modular data. 

Proof. For (1), see [BKj Eq. (3.1.2)] on page 47. For (2), it is [EH Theorem 
3.1.19] found on page 57. Formula (3) from |Baj for RCFTs can be generalized to 
MTCs. (4) follows from Theorem 12.51 □ 

Proposition 12.31 (2) implies that the 0j are actually roots of unity of finite order, 
which is often referred to as Vafa's Theorem. But from example I2.2[ we know 
that this is not true for general modular symbols, in particular Q(8i) might not 
be algebraic for modular symbols. This leads to: 

Definition 2.4. Given a modular data (N;S,T), let %n be the Galois field 
Q(sij,D,9i),i,j G £. Then the Galois group of %n over Q will be called the 
Galois group of the modular data. 

Theorem 2.5. (1) (de Boer-Goeree theorem): The Galois group of a modular 
fusion rule is abelian. 
(2) The Galois group of a modular data is abelian. 

By the Kronecker- Weber theorem, there is an integer m such that %n C Q(Cm)j 
where ( m = e^r '. The smallest such m for %n is called the conductor of %n, and 
the order of T always divides iV (we intentionally build iV into the notation %n)- 
The Galois group of Q(Cjv) is the cyclic group of units I such that gcd(Z, N) = 1. 
Each I acts on % N as the Frobenius map 07 : — > ( l N . Consequently, 07 (T) = T l 
and &i(S) = SP a , where the signed permutation matrix P a corresponds to the 
Galois element a in the Galois group of the modular fusion rule. 

It is known that the fusion algebra of a rank=n MTC is isomorphic to the func- 
tion algebra of n points. A Galois group element a of the associated modular 
fusion rule induces an isomorphism of the fusion algebra. It follows that a deter- 
mines a permutation of the label set. When we have only a modular fusion rule, 
the two algebra structures on the fusion algebra a priori might not be isomorphic 
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to each other. But still a Galois group element of the modular fusion determines 
a permutation of the label set and the de Boer-Goeree theorem holds. Actually 
what we are using in this paper are identities among modular S entries up to some 
parity signs e» )<T = ±1 associated to each Galois element a. Such parity signs first 
appeared in |CGj for Galois automorphisms of Q(\ij,D). 

First we note the following easy, but very useful fact that the ordered set of 
eigenvalues of iVj determines the label i uniquely. 

Proposition 2.6. There do not exist indices j ^ k such that Ajj = \^ for all % 
for any modular fusion rule (N; S) . 

Proof. If there were such indices, then the dot product of rows j and k of S would 
be D 2 = J21=o \$i,j\ 2 > 0, a contradiction. 

□ 

Except (5), the following theorem is contained in [CGJ. 

Theorem 2.7. Let G be the Galois group of a rank=n modular fusion rule (N;S). 
Then 

(1) the simultaneous action of the Galois group G on the set A = { A^ } gives 
an infective group homomorphism i : G — > S n , where S n is the permutation 
group of n letters; for a G G, i(a)(i) is the associated element in S n ; 

(2) for any a G G, the matrix P a = c/ cr (o)5' _1 cr(5') is a signed permutation 
matrix; furthermore, the map o — > P a gives a group homomorphism from 
G to the signed permutation matrices modulo ±1 which lifts i; 

(3) for each a G G, there are 6{ ja = ±1 such that 

(2-4) a (Sj,k) = ^a(k),aSj,(T<k)- 

«<t(0) 

Moreover, 

(2.5) Sj t k = e<T(jV e fc,<^OV -1 (fc)> 
and 

(2.6) £<T- 1 (k),<r- 1 — e<T(o),0-eo,crefc, CT ; 

(4) the Galois group G is abelian; 

(5) If n is even, then 11™=$ ej. CT = (— 1)°". If n is odd, then D G K, and 
a(D) = e a ■ ^ where e a = ±l. We have U?=o U, a = ^ • (~l) a ■ 

We are going to use a for both the element of the Galois group G and its 
associated element of S n . When o G G applies to a matrix, a applies entry-wise. 

Proof. Let K = Q[{^i,j}o<i,j<n-i] be the Galois extension of Q generated by the 
eigenvalues of all the A^j and let G be the associated Galois group as above. The 
action of G on the eigenvalues gives an injection G ^ S n x S n x ■ ■ ■ x S n , where 
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there are n — 1 factors. Note that we have not assumed the Ni have distinct 
eigenvalues, therefore this map is not necessarily unique and is not necessarily a 
group homomorphism. This is not a problem as we will resolve the ambiguity 
shortly. Just fix one such map for now. Let (a%, cr 2 , • • • , u n -i) denote the image 
of cr g G under this injection. Note that a priori, there is no relationship between 
the <7j. Let A, be the diagonal matrix with diagonal entries Ay, so 

Ni = SAiS' 1 . 

Let P ai = (5i=ai(j))a<i,j<n-i be the permutation matrix corresponding to cr^. Since 
a (^i,j) = ^wO')' we nave c"(Ai) = P~^kiP ai . Since is rational we have 

SA.S- 1 = N = a(N t ) = aCS)P; i l KPaA~SY 1 - 
Rewriting this gives 

A^-V^P- 1 ] = [5-V(5)P- 1 ]A i , 

Hence B it(7 = S^ 1 cr(S)P^ i i commutes with Aj. It follows that B i)(J is block diagonal, 
with blocks corresponding to the equal eigenvalues of iVj. In formulas, if the (J, k) 
entry of B ia is nonzero, then Ay = Let S~ x o{$) = B i a P a . = C a . Note two 
facts, if the (J, k) entry of C a is nonzero, then the (j, (Ji{k)) entry of B i)(J is nonzero 
and hence X t j = Xi iai (k)- The second fact is that C a (as the notation suggests) does 
not depend on i, only on a. Suppose C a has 2 nonzero entries in column k, say 
the (j, k) and (£, k) entries. Then Ay = Aj jCTi (fc) = \ t i for all i, contradicting the 
Proposition 12.61 above. If a row or column of C a is all zeroes, then det(C cr ) = 0, a 
contradiction. Hence C a has exactly one nonzero entry in every row and in every 
column. Thus there is a unique permutation a £ S n and a diagonal matrix B a 
such that C a = B a P G . Note that we are now using a for both the element of the 
Galois group and its associated element of S n . Note that 

C aa , = S^aa'iS) = 5-V(5)«7(5 _ V(5)) = C a a{G a ,). 

From which it follows that the map G — ► S n is a group homomorphism. Thus we 
have proved that the simultaneous action of the Galois group G on the eigenvalues 
Aj.-, of Ni for all i gives an injective group homomorphism G — ► 5 n . 

Note that the squared length of column zero of >S is D 2 = ^7=0 d% which must 
be equal to the squared length of column cr(0). Hence 

n—l /n—1 \ 

D 2 = J2 rf a(0)\>(0) = ^(0)°" Z! A ?0 = d l(0) a ( D2 )- 

i=o V i=o / 

Rewriting gives 




It follows that G acts in the same way on the quantities {dj/D 2 }. The Verlinde 
formulas (12.31) encode the symmetry of the Ni matrices, and give us the complete 
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symmetry under interchanging the last n — 1 N{ and simultaneously reordering the 
last n — 1 rows and columns of all matrices. Thus „■ is invariant under G and 
hence is necessarily rational if we define it first to be only in M. 

Transposing the identity S' _1 o"(S') = C a and inverting this identity gives the two 
equations = Gj and 



D 



2 



T 



c- 1 = aisy's = -^a&s- 1 = 4 (0) g, 

Hence the matrices d a ^C a and d a ^B a are orthogonal. Since B a is diagonal it 
follows that 

/e 0>CT ■•• \ 
e liCT ••• 



5, 



^(x(0) 



V Cn-W 

for some choices of ej )tT = ±1. The map a i— > d^Go- gives a group homomorphism 
from G to the signed permutation matrices modulo ±1 which lifts the homomor- 
phism i of (1). 

Rewrite the definition of C a as u(S) = SB a P a . Picking out the (j, k) en- 
try, we have cr{sj^) = j^£a(k),oSj,a(k)- Moreover, since the left hand side is 

symmetric we get SB a P a = P a : 1 B a S. In coordinates this condition becomes 
Sj,k = ^k^a(j)S a -(j),cT- 1 (k)- Consider the action of G on pairs (j, k) defined by 
a x (j, k) i — > (cr(j), cr _1 (A;)). Then we see that \Sj t k\ is constant on orbits of this 
action. 

To see identity (12. 6p . we apply cr -1 to identity (12.41) and compare with identity 
(J23]). Note that s a -i {0):(7{0) = e a{Q ),a^o,a by identity (J53D. 

Given a x ,a 2 E G, consider first <7 2 0i(%,fc) = ^{j^e^^Sj^k)) 

Then consider 

0"lO"2(Sj,A ; ) = Oia 2 {Sk,j) = a l{~^ e a-2(i),a- 2 5fc,a- 2 0)) = 

"<T 2 (0) 

{ 1 ~ 1 

°H , e (T2(j),(T2 S V2(j),k) , x e CT20'),o-2 e cri(A:),(Ti S o-2(i),o-i(fc)- 

«CT2(0) «fTl(0) A O-2(0),(Tl(0) 

Hence cricr 2 = cr 2 eri using diXjj = dj\ it j, i.e. G is abelian. 

Suppose now that the rank n = 2r is even. Then det(S') 2 = D 2n hence 
det(S') = ±D 2r . Since the determinant is a polynomial in the entries of the matrix 
det(cr(5')) = ±<j(D 2 ) r , with the same sign as det(S'). Hence det(S'~ 1 cr( 1 S')) = d~JL. 
Since det(G CT ) = d~Zs(—l) a Uj=o we conclude UjZo e j,<? = For °dd rank 

n = 2r + 1, det(S) = ±D 2r+1 , hence D e K. Hence a(D) = e a D/d a(0) ,e a = ±1 
and one gets the formula IljZo e j,<? — e o-( — 
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□ 

Note that the resulting equations ( 12. 5ft for the entries §j t k are unchanged if we 
replace B a with —B a . We will use this to assume eo = 1 below. 

Next we will use the fact that the 9i e U(l) to produce a series of twist inequal- 
ities on the entries of S. 

Theorem 2.8. Given a modular symbol (CNT; S,T) and S is a real matrix, then 

(1) 2maXj s| • < D\Sjj\ + D 2 for any j . 

(2) If j ^ k, then D<^ J27=o \hAk\- 

Proof. Rewrite the twist equation as TSTST = D + S. Then taking the (j, k) entry 
of this formula gives 

n-1 
i=0 

Since \D + \ = D and \6i\ = 1, the largest of the n + 1 numbers |si,jSi,&|, < 
i < n — 1, and -D|Sj.fc| must be at most the sum of the other n. If j = k, then 
J2i$ij = D 2 > D\Sjj\. Hence this inequality is trivial unless the largest is one of 
the first n and we get 

n-1 



2maxs 1j < D\hi\ + £ ~ s lr 



i=0 



If j ^ k then X)j Si,jS-i,k — and the nontrivial case is 



i n—1 



-° ^ p — r Yl \hihk\- 



\ s j,k 



i=0 



We will refer to these as the twist inequalities. 

Suppose er G G corresponds to signs Cj as above. We drop a for notational 
easiness. Multiply the identity above by £a(j)/(QjOati)), set k = cr(j), and sum over 
j. The result is 

n—1 n—1 n—1 ~ 

Z Mi) Z Qihihati) = D + J2 G a a h ° 3 • 

j=0 i=0 j=0 °j°<?U) 

Interchanging the sums and using the fact that Si )(7 (.j) = £a(j)Ca(i)S<j(i),j gives 

n — 1 n—1 n—1 ~ 

i=0 j=0 j=0 a 3 u <r(j) 
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By orthogonality of the rows of S, the innermost sum on the left is zero if % ^ o~(i) 
and D 2 = D + D_ if i = a(i). Hence 

n—l ~ 

j=0 u j u °{i) i; a (i)=i 

If a is fixed point free, then £?-n ^dpfte® = 0. 

□ 

3. RANK=2 AND 3 MODULAR S MATRICES 

In this section, we determine all possible modular S matrices for rank=2 and 3 
modular symbols. The rank=3 case first appeared in |CP] , but our proof is new. 

Theorem 3.1. The only possible rank=2 modular S matrices of some modular 
symbols are 



(1) 



where e 2 

(2) 



where tp 2 = 1 + (p. 



1 e ' 



'1 ^ 



Proof. Since all labels are self-dual, S is a symmetric real unitary matrix of the 

form \ \ . The fusion matrix Ni is of the form ^ , so we have d 2 = 

1 + md. Simplifying D + D_ = D 2 leads to 9 + 9~ L = 1 — d 2 = —m ■ d. Since 
9 E U(l), so \md\ < 2. If d > 0, then d = m+v ^ ? + i , hence m = 0, 1. If d < 0, 

then d = m ~ v ^ ?Ti , hence [Q(0 + : Q] < 2. It follows that 9 = eT for some 
(PiO) = 1; an d 9 is one of {1,2,3,4,5,6}. Direct computation shows there are 
no integral solutions p,q for 2cos(^ E ) = — m ■ m ~ v ^" 2 + 4 except for q = 2, 5 and 
m = 0, 1. 

□ 

Theorem 3.2. T/ien i/ie onZy possible rank=3 modular S matrices of some mod- 
ular symbols up to permutations are 

(1) 

/l e e 

,2 



6 W W 
t UJ 2 UJ 
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(2) 



where e 2 = 1, and u> 3 = l,u> ^ 1. 



/I 

d 



V 



where d 2 = 2. 



(3) 



d 









-d 




d\ 




-d 2 




1 


— d\ 



fl 

d\ 
\d 2 

where d\ is a real root of x 3 — 2x 2 — x + 1 and di = di/{di — 1) which is 
a root of x 3 — x 2 — 2x + 1 . T7ie largest di 
and d 2 = 2 cos(tt/7) = 1.801937736 .... 



2cos(tt/7) 
2cos(tt/7)-1 



2.246979604 



Proof. The case of non-self dual case is given in Appendix |A] Hence we assume all 
fusion rules are self-dual, so S is a real, symmetric, unitary matrix up to the scalar 
D. It follows that the fusion matrices iVj's are commutative, symmetric, integral 
matrices. One approach to proving the theorem is to analyze case by case for the 
Galois groups of fusion rules G = 1, Z 2 , Z 3 . This strategy will be fully exploited 
in the rank=4 case in the next section. Instead we will argue directly from the 
^-matrix in this section. 

The fusion matrices Nx,N% are symmetric, and NiN 2 = N 2 Ni. Therefore, they 
can be written as 

/0 1 N 
Ni = 1 m k 
\0 k I 

and 




No 



such that 

l + ml + kn = k 2 + I 2 . 
There characteristic polynomials are 

Pl (x) = x 3 -(£ + m)x 2 + {mi - k 2 - l)x + t = 

and 

p 2 (x) = x 3 - (k + n)x 2 + {nk - I 2 - l)x + k = 0, 

respectively. 
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Next we turn to the S matrix, which is of the following form: 

/ 1 di d 2 \ 

S = d\ Su §12 

\d 2 §12 §22 / 

Orthogonality of the columns of the S matrix translates into the equations 

di + di§u + d 2 s 12 = 0, 
d 2 + c?i§i2 + d 2 s 22 = 0, 
did 2 + §i2(§n + §22) = 0. 

The first two equations give §n = —1 — ^2 §12/^1 and §22 = — 1 —d\ §12 /d 2 - Plugging 
these into the third equation gives 

(d\ + ^2)^12 + 2did 2 s 12 — d\d\ = 

hence 

d±d 2 

S12 



1±D 

Thus 



1±D' 

and 

dj 



hi = -1 - 
§22 = -1 



dl 



1±D 

Thus the eigenvalues of Ni are 

§12 di §11 1 d\ 

di, b — — = — — — , and c — 



d 2 1±D' d 1 d 1 d 1 (l±D) 

and the eigenvalues of N 2 are 

, §12 d 2 , , 1 d\ 

«2, e = — = — — — , and / = — — 



di 1±D' J d 2 d 2 (l±D) 
We compute 

dib + d 2 f = die + d 2 e = be + ef = — 1. 
Since dibc = — £ and d 2 ef = —k, these are equivalent to 

£ /c £ /c £ fa 

c e b j di GS2 

Also note that 

die = d 2 b. 

Let's deal with the case where £ = first. Then we have k 2 = kn + 1. Hence 
k — 1 and n = 0. Thus the eigenvalues of A^ 2 are 1,1, and -1 and the eigenvalues of 
Ni are (m+y/m 2 + 8)/2, 0, and (m — ^/m 2 + 8)/2. Since Aq has eigenvalues c?i, b, c, 



16 



ERIC ROWELL, RICHARD STONG, AND ZHENGHAN WANG 



and d\ ^ 0, hence c = which implies m — 0. This gives (k, £, m, n) = (1, 0, 0, 0) 
and 

(\ d 1 
S= d 
\1 -d 

where d 2 = 2. The case k = gives essentially the same solution, so we will 
henceforth assume £ and k are positive. Since pi(^) = — k 2 £ < and pi(0) = £ > 0, 
we see that the largest root of p\ is > £, one of the remaining roots is in (0, £) and 
the other root is negative. Similarly the largest root of p 2 is > k and the other 
roots are in (0, k) and (— oo, 0). 

Case 1. The polynomial pi(x) is reducible. 

Since d\ > £, d\ cannot be an integer. Thus p\ must split into a linear and 
an irreducible quadratic. Thus Q[di,D] is a quadratic extension of Q. Hence 
Q[di, d 2 , D] has degree 2 or 4 over Q. Thus p 2 is also reducible and also splits as 
into a linear and an irreducible quadratic. Since the £/b + kj f — £/c + k/e = 1, 
the integral roots must be either b and / or c and e. Without loss, we may assume 
the integer roots are b and /. Let 

d 2 = a + (5\fs and e = a — (3y/s 

for rational (in fact integer or half-integer) a and (3 and integer s. Then since 
die = d 2 b and c is the conjugate of di we have 

di = b — F and c = 6 — =. 

ct — pv s a + Pv s 

Hence £ = -d x bc = -6 3 . Since / = -k/(d 2 e) = -k/(a 2 - (3 2 s) and -b 2 = £/b = 
1 — k/ f = a 2 + 1 — j3 2 s. Therefore solving 1 = k/d 2 + £/d 1 for k gives 

k = d 2 - £^- = a + + b 3 l = a + (3y/s - (a 2 + 1 - (3 2 s) (a - 0y/s) 
d\ b 

= -a(a 2 - (3 2 s) + (3[a 2 + 2 - (3 2 s]Vs. 

Since k is an integer, this forces a 2 — f3 2 s = —2, hence b 2 = 1. Since £ > 0, 
this means £ = 1 and 6 = —1. Also from the equations above we get k = 2a, 
d 2 = a + \Jot 2 + 2, e = a — V" 2 + 2, / = a, d 1 = a 2 + 1 + ay/ a 2 + 2 = d^/2, 
c = a 2 + 1 - aVa 2 + 2, and £> = a 2 + 2 + «V« 2 + 2. Thus 

= x 3 - (2a 2 + l)a; 2 - (2a 2 + l)x + 1 

and 

p 2 (x) = x 3 - 3ax 2 + (2a 2 - 2)x + 2a. 
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Thus (k,£,m,n) = {2a, 1, 2a 2 , a) and 



/ 1 a 2 + 1 + aVa 2 + 2 a + V« 2 + 2 



2 



1 + a\Ja 2 + 2 1 -a-Vo? 



\ a + y/a 2 + 2 -a - \/a 2 + 2 a 2 + aV« 2 + 2, 

Note that n = a must be a non-negative integer. Setting a = gives the example 
found above again. Thus we may assume a > 1. Since d\ = d 2 /2, the equation 
for the #'s is 

|1 + 9 2 d\ + (1/4)^4| = D = 1 + (1/2)4 
If a solution did exist, then we would have 1 + (l/2)d\ > (l/4)rf| — d\ — 1, hence 



17 > (d 2 . - 3) 2 or d 2 < V3 + = 2.66891 .... Since d 2 > 1 + V% = 2.73205 
this cannot occur. Thus these S matrices, for positive a, do not give a modular 
symbol. 

Case 2. The polynomial Pi(x) is an irreducible cubic. 

By Case 1, we see that p 2 (x) is also an irreducible cubic. Then there must be 
a Galois symmetry o with a(di) = b, a(b) = c and a(c) = d±. Hence cr(d 2 ) = f, 
cr(f) = e, and a(e) = d 2 since these roots of p 2 pair with the corresponding roots 
of p\. Applying a to the identity d\t = d 2 b, gives d 2 b = fc. Thus we must have 
fc = d l d 2 /(l±D). Since 

d x d 2 d 2 + d 2 1 d x d 2 D 2 ±D 

J c = m i m9 + j j /-, i + 3~T = m i m9 + 



(1±£>) 2 d l( i 2 (l ±D) did 2 (1±L>) 2 ^(liD)' 
we compute 

1±D= K J fc= 1± 



and hence 



Since D > 1, we get that 



did 2 " d\d\ 



did 



2 



1. 



d\d 2 

Sl2 = = ±1 

12 1±D 

and hence 6 = ±l/d 2 and e = ±l/di. Thus d\ and rf 2 are units in the ring of 
algebraic integers. Hence k — £ — 1 and hence m + n = 1. Without loss we may 
assume m = 1 and n = 0. Then pi(:r) = x 3 — 2x 2 — re + 1 and p 2 (x) = x 3 p\(l/x) = 
x 3 — x 2 — 2x + 1. Then one computes 



18 



ERIC ROWELL, RICHARD STONG, AND ZHENGHAN WANG 



b = 1 — 1/di, C 

f — 1 — di and 



-l/(di - 1), d 2 = d 1 /(d 1 - 1) = 1.801937736 . . . , e = l/d u and 



S 




□ 



4. RANK=4 MODULAR S MATRICES 

First we introduce the following notation. For an integer m, define 

m + y/m 2 + 4 

0m = g ' 

that is, m is the unique positive root of x 2 — mx — 1=0. Note that any algebraic 
number whose only conjugate is —1/0 must be m for some integer m. Also 
note the only rational m is O = 1. 

Theorem 4.1. TTie on/y possible rank=4 modular S matrices of unitary modular 
symbols up to permutations are 

(1) 



(2) 



where u = ±i; 



(3) 



(4) 



A 


i 


i 


1\ 


i 


i 


-i 


-1 


i 


-i 


a; 




V 


-i 


u; 


u J 


n 


i 


1 


1 \ 


i 


i 


-1 


-1 


i 


-i 


1 


-1 


V 


-i 


-1 


1/ 


n 


i 


1 


1 \ 


i 


-i 


1 


-1 


i 


i 


-1 


-1 


V 


-i 


-1 


1 / 


n 




1 






-i 




-1 


i 




-1 




VP 


-i 




1 / 



where (p = ^y^- is the golden ratio; 
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(5) 



(6) 





( 1 <P 


<p <p 2 \ 








up -1 


v 2 -<p 








<p V 2 


-1 -if 


J 








-if 1 ) 






1 


d 2 -l 


d + 1 


d 




d 2 - 1 





-d 2 + l 


d 2 - 


1 


d+1 


-d 2 + l 


d 


-1 




d 


d 2 - 1 


-1 


-d- 


V 



where d is the largest real root of x 3 — 3x — 1 . 

Proof. The non-self dual case is treated in Appendix [A], so we will assume that S 
is real in the following. Since the fusion coefficients n^j are totally symmetric in 
i,j and k for self-dual categories, we will instead write in what follows. For 
notational easiness, when the Galois group element a is clear from the context, we 
simply write ej )(T as Identities for e« and Sjk that are not referenced are all from 
Theorem 12.71 All the twist inequalities are from Theorem 12.81 



Case 1. G contains a 4- cycle. 



By symmetry we may assume a = 

€k£a(j)Sa(j),<r-Hk) ancl ^1^3 = ~l give 

/I 

3 = 



(0 1 2 3) G G. The conditions Sj^ 



\d 3 



d\ 
e 1 e 2 d 2 
-e 2 d 3 



\ 



d 3 

e 2 d 1 
-eie 2 d 2 J 
we may assume e 2 



d 2 
-e 2 d 3 
-1 
e 2 di 

By symmetry under interchanging Aq and iV 3 , we may assume e 2 = +1. Note 
that o~ 2 (d 2 ) = X 2t2 = —l/d 2 ^ d 2 . Hence the characteristic polynomial p 2 of N 2 is 
irreducible. Since cx 2 (di) = —d 3 /d 2 < 0. Hence cr 2 (<i 1 ) ^ d\. Thus p\ is irreducible. 
Since ei/d 3 is a root of p±, it follows that p 3 is also irreducible. 

We see that Ai,i = t\d 2 jd\^ \\ )2 = —d 3 /d 2 , and Ai,3 = e\/d 3 . In particular 
diAi^Ai^Ai^ = —1. Orthogonality of the rows of S is equivalent to 

di + e\d\d 2 — d 2 d 3 + e\d 3 = 0, or 
1 1 



+ Ai 3 — Ai i + 



d\ Ai f 2 
Write pi(x) = x 4 — Cix 3 + c 2 x 2 + c 3 x — 1. Thenp 4 (x) = x 4 — eic 3 x 3 — c 2 x 2 + eiC\x — 1. 
Note that c\ = Trace (Aq) > and e x c 3 = Trace(A 3 ) > 0. Multiplying together 
the orthogonality condition above and five of its formal conjugates gives 

128 + (c 2 3 - c 2 ) 2 - 16cic 3 + 12(c 2 - c\)c 2 + 32c 2 2 = 0. 
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This equation forces c\ and C3 to be even. Let A = c\ — C3 and E = c\ + C3 (hence 
E and A are even and congruent mod 4). Then solving the quadratic equation 
above for c 2 we see that we must have (E 2 — 32) (A 2 + 32) to be a square and 

c 2 = 3AS=fc V( s i6 32 ' )( ' A 't M, it follows that |E| > 6. If E and A are multiples of 4, 
then we see they are multiples of 8 and either sign gives an integral c 2 . If E and A 
are both 2 mod 4, then there is a unique choice of the sign for which c 2 is integral. 

The Galois group of p 1 must be Z/4Z, otherwise it would contain the (0 1)(2 3). 
Applying this to the orthogonality identity above gives 1/Ai, 3 + d\ — Ai, 2 + 1/Ai,i- 
Multiplying this by the original identity gives - + diXi^ = Ai 2 + Ai,iAi )2 . 

Hence d±Xi^ = (Ai,iAi, 2 ) ±:L , either of which contradicts the product of all four roots 
being —1. In particular, p 1 cannot have complex roots, since complex conjugation 
would give a transposition in the Galois group. Applying a to the orthogonality 
identity gives d\ + 1/Ai,i = Ai, 2 + \ jX\^. 

We know from the preliminary discussion that all three of the resulting N{ 
matrices will be rational. Define P = 16 ^ 2 2 ^ 3 3 AS = ±^ ^+32 "' then we com P u te 

5ci - 3c 3 ci - c 3 
= — g g— P, 

"1,1,2 = ei(P - 1), 

/C1 + C3 C1-C3 
"1,1,3 = ei — — P 



Cl + c 3 Ci - c 3 
"1,2,2 = + -^P' 

"1,2,3 = P 

Cl + C 3 Ci - c 3 
"1,3,3 - g g F, 

n 2 , 2 , 2 = ei (2c 2 - - 2p) , 

( Cl + C3 , C1-C3 p \ 

"2,3,3 = ei(P + 1), and 

/5c 3 -3ci C1-C3 
"3,3,3 = ,i g— P 

Recall that the riij,fc must be nonnegative integers. This restricts the q. First 
looking at ni, 2 ,3, we see that P must be a positive integer. Hence c 2 must be given 
by the upper sign. This condition in fact guarantees integrality of all the n^j^- 
(The additional factors of 2 in the denominator cancel out if c 2 is integral and can 
be ignored.) Integrality of P severely restricts A, since it requires all odd prime 
factors of A 2 + 32 to be congruent to 1 mod 8 (since 2 and —2 are both squares 
mod any such prime). In particular either A = or |A| > 6. Since P ^ 0, we see 
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that £ 2 > A 2 + 64, hence |S| > |A|. Thus c 3 must be positive and E > 0. Since 
we saw above eic 3 > 0, we see that e x = +1. Thus rewriting the orthogonality 
relation gives d\jd 3 = (d 2 — l)/(d 2 + 1). 

The twist inequality coming from the (0, 3) entry reads 

D <(i + *) ( i +<w . 

Plugging in the preceding identity, simplifies this to D < 2d 2 . Rearranging gives 
3c? 2 > d\ + d\ + 1 > d\ + d\ and plugging in the identity d\ = d 3 (d 2 — l)/(d 2 + 1) 
yields 

/ d 3 V U\ 3 



< — 2 - <^ 



\d 2 + lj 2(d| + 1) 2 

To see why this is helpful, expand the equations Trace(iVj) = q for i — 1,3 and 
use the identity above to eliminate d\. The result is 

d\ - 2d 2 - 1 , d| + 2d 2 - 1 1 
Cl = , 7 , . 7T rf 3 + "j : X' alld 

^ + 2d 2 - 1 , d| - 2d 2 - 1 1 

c 3 = m 1 -1 \ 3 



d 2 (d 2 + 1) d 2 — 1 d 3 

Subtracting these gives 

, d 3 d 2 + l 

A = ci - c 3 = -4-— + 2——. 

a 2 + 1 0203 

Hence A > -4^_ > -4^3/2 > -4.9. However, we saw above that either A = 
or |A| > 6. It follows that A > 0. Since £ > 6 and c 2 > 3AE/16, it follows that 
c 2 > A = ci — c 3 . Thus p 3 (l) = Ci — c 2 — c 3 < 0. Thus rf 3 > I. Hence we see 
A < 2{d 2 + \)/{d 2 d 3 ) < 4. It follows that A = 0, i.e., c x = c 3 . 
Since A = 0, £ = 2ci is a multiple of 8 and 

^'- 2 f£iV = -i. 



V32/ 

In this case the characteristic polynomials become 

Pl (a;) = a; 4 - cix 3 + 2^2(a/4) 2 - lx 2 + c x x - 1, 



P2( x ) = X A - 4J2(ci/4) 2 - 1^ - Qx 2 + 4J2(ci/4) 2 - lx + 1, and 



p 3 ( x ) = X A - Cl a; 3 - 2 v /2(ci/4) 2 - lx 2 + cix - 1. 
In particular pi(x) > for x > c\. Hence d\ < c\. We have 

p 2 (x) = (x 2 - tix - l)(x 2 - t 2 x - 1), 
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where ti > > t 2 are the two roots of t 2 - 4^/2(ci/4) 2 - It - 4 = 0. Since the 
larger root of x 2 — tx — 1 is an increasing function of t, d 2 must correspond to t\. 
Hence 

d 2 = h + J- > h = 4 v /2(c 1 /4) 2 - 1 + i > 4 v /2(c 1 /4) 2 -l. 

«2 H 

In particular, d 2 > 4 since the square root above is integral. Finally the twist 
inequality coming from the (0, 2) entry reads 

d 1 d 3 \ 



D < 2 1 + 



do 



Squaring and using the identity d 3 = di(d 2 + 1) / (d 2 — 1) to eliminate d 3 gives the 
inequality 

4(rf 2 + l) 2 d{ - 2d 2 {d\ - Ad\ + d 2 + A)d\ - d\{d 2 - l) 2 (d\ - 3) > 0. 
Dividing through by A(d 2 + l) 2 d\ and rearranging gives 

The right hand side of this inequality is an increasing function of d 2 for d 2 > 4 
and a decreasing function of d\, hence we may replace d\ by its upper bound C\ 
and d 2 by the lower bound above. The result is 



4560 - 2138c 2 + 264c^ - 8c? + (608 - 276c 2 + 32^)^? - 16 > 0. 

Since the coefficient of the square root is nonnegative it follows that 

4560 - 2138c 2 + 264c? - %c\ + (608 - 276c 2 + 32c?)c 1 v / 2 > 0. 

This polynomial in c\ is negative for c\ > 6\/2, hence c\ < 8. The only multiple of 
4 in the range 3 < c\ < 8 for which the Pell equation above is satisfiable is c\ = 4. 
This gives c\ = c 3 = 4, c 2 = 2, P = 32. However plugging in shows that the twist 
inequality D < 2(1 + did 3 /d 2 ) does not actually hold in this case. Thus there are 
no solutions in this case. 



Case 2. G is the Klein 4-group. 

Let <7i, 02, and 03 be the elements of G which correspond to (0 1)(2 3), 
(0 2)(1 3), and (0 3)(1 2), respectively. Let B ai correspond to signs e { with 
£it2t3 = 1 and B a2 correspond to signs 5i with ^i^^a = 1- Then using the usual 
identities gives 



/I 

d\ 
d 2 
\d 3 



di 
ei 

6i6 2 d 2 



d 2 
e 2 d 3 

$2,2 
§2,3 



d 3 \ 

6i6 2 d 2 

§2,3 
ei§2,2/ 



/I 

d\ 
d 2 
\d 3 



d\ 
§1,1 
$id 3 

§1,3 



d 2 
$id 3 
S 2 
S\5 2 di 



d 3 \ 

S l,3 

5i5 2 di 
S 2 h,iJ 
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Comparing these we see e 2 = 81, §1 



e\, and 82,2 = 82, hence 



/l di 

d\ ei 

\d 3 e 1 e 2 d 2 



d 2 
e 2 d 3 

82 
e 2 8 2 di 



d 3 \ 
eie 2 d 2 
e 2 8 2 d\ 
£182 J 



Orthogonality of the rows of S gives the three conditions 

(1 + €i)(di + e 2 d 2 d 3 ) = (1 + 8 2 )(d 2 + e 2 d 1 d 3 ) = (1 + €i8 2 )(d 3 + eie 2 did 2 ) = 0. 

Suppose ei = +1, then we see d\ = —e 2 d 2 d 3 , hence e 2 = —1 and the remaining 
orthogonality relations become d 2 (l + 8 2 )(l — d 3 ) = d 3 (l + 8 2 )(l — d 2 ). We cannot 
have d 2 — d 3 — 1, since this would make d± — 1, hence 8 2 = — 1. This gives 

/ 1 rf 2 C?3 c?2 d 3 ^ 

^2^3 1 — d 3 —d 2 

d2 -d 3 -1 d 2 d 3 

\ d 3 —d 2 d 2 d 3 -1 J 

The eigenvalues of N 2 are rf 2 and —l/d 2 each with multiplicity 2. Hence d 2 = 4> m 
for some integer m. The eigenvalues of N 3 are c?3 and —l/d 3 each with multiplicity 
2, hence d 3 = (f> n for some integer n. So 

/ 1 m n m n \ 

-</>« -1 



5 = 



'mrn 



V 0n -0m 0m 0n ~1 / 

The resulting iVj matrices are necessarily rational, but in this case they are all 
integral, namely, 

/o 1 0\ 

1 mil m n 



JVi = ^2^3 = iV 3 iV 2 







m 
n 




1 



1 



/o 

m 

1 
VO 1 



1 



1 


m 

0; 



and A^3 



/0 l\ 

n 1 

10 

VI n) 



Note that nonnegativity of the entries forces m, n > and hence m , 0„ > 1. The 
strongest twist inequalities are the (0, 1) and (2, 3) cases which give D < 4 or 
(0^ + 1)(0^ + 1) < 16. This gives, up to symmetry, the solutions (m,n) = (0,0), 
(0, 1), (0,2), or (1, 1). These are all excluded since the resulting Galois group G 
is at most Z/2Z. (These examples will return when we look at smaller Galois 
groups.) 
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Case 3. G contains a 3-cycle. 

Since we can exclude cases 1 and 2 above, the image of G in S4 cannot be 
transitive. It follows that G must fix the point j not on the 3-cycle. Thus Ajj is 
rational (hence integral) for every %. Up to symmetry there are two cases for the 
3-cycle. We could have a = (1 2 3) or a = (0 1 2). If a = (1 2 3), then the d { 
are integral. The identities Sj^ = taij^kSaij)^- 1 ^) and eie 2 e 3 = 1 give = 1 for 
all % (since di = s 0:i = ejS ,j+i — Q^i+i for 1 < i < 2) and 



Orthogonality of the columns of S gives §1,1 + §2,2 + 53,3 = — 1 and Si,is 2)2 +s 2 , 2 s 3 , 3 + 
§3,351,1 — —d\. The first of these gives — \jd\ = A^i + A 12 + Ai )3 from which we 
see —1/di is an algebraic integer. Hence d\ — 1 and A^j < 1. The second equation 
gives Ai ; iAi i2 + Ai )2 Ai )3 + Ai^Ai^ = —1. Hence Ai^, Ai j2 , and Ai j3 are the three 
roots of g{x) = x 3 + x 2 — x + n for some integer n. This cubic must be irreducible 
and have three real roots all less than 1. Irreducibility excludes n = and n = — 1. 
For the roots of g to be less than 1, we must have g(l)>0orn+l>0. Hence 
n > 1. However, this results in complex roots. Thus this case gives no solutions. 



Thus we must have a = (0 1 2) and A ij3 is integral for all i. The identities for 
Sj,h give 



Since a(d 3 ) = A 3j i = e 2 d 3 jd x and <J 2 (d 3 ) = €ie 2 d 3 /d 2 , we must have a(d 3 ) 7^ d 3 . 
(Otherwise e\ = e 2 = d\ = d 2 = 1 which fails.) Thus d 3 is a root of an irreducible 
cubic g(x) = x 3 — C\x 2 + c 2 x — c 3 and e 2 d\ and t\t 2 d 2 are ratios of roots of g. 
If g had Galois group S3, then the ratios of the roots of g would be roots of an 
irreducible sextic. Thus g has Galois group Z/3Z and G = {l,a,a 2 }. Note that 



S 



/ 1 di di d 1 \ 

di §1,1 §3,3 §2,2 

d\ §3,3 §2,2 Si,l 

\«1 §2,2 §1,1 §3,3/ 



5 



/ 1 C?x d 2 ^3 \ 

d 1 t x t 2 d 2 ei e 2 (i 3 

c?2 ei £2^1 ^i^d 3 

\d 3 e 2 d 3 tie 2 d 3 s 3i3 / 



Cl = d 3 + A 3j i + A 3i2 



d x d 2 



[d\d 2 + e 2 c?2 + eie 2 di), 



c 2 — 6^3,1 + A 3) iA3 i2 + A 3j2 0? 3 



(e 2 o?2 + ei + eie 2 di), and 



C3 — ^3A3,lA 3)2 — €i 
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Orthogonality of the columns of S gives 

Cl 62^1 + ^ie 2 d 2 + t\d\d 2 

— = ^ = — 1, and 

c 2 1 + e 2 di + e x e 2 d 2 s 3 , 3 x „ 

— — - — — — — — a 3i 3 e ^. 

c 3 d 3 a 3 

Thus = — ca; 2 + ncx + c for integer n, c. Since g has Galois group Z/3Z, 



1 



5 2 = — discr(^) = (n 1 + 4)c 2 - 2n(2n 2 + 9)c - 27 



must be a square. The resulting n^j^ are 

n 1 , 1 ,2 = ^^yM + nc-2n 2 + 3), 
721,1,3 = 2(^+3) + ^ + + 

ni,2,2 = 2RT^) (5 + nc - 2n2 + 3) ' 

ni,2,3 = ^h (3n - c) ' 

n 1A3 =2j^^(5-nc + 2n 2 -3), 
»2A2 = — (5 + nc + l) + -^, 

"2,2,3 = 2( - n 2 + 3) ( n<5 + + 2 ) C + 3n )> 

n 2 s 3 = , 1 2 — r(— 5 — nc + 2n 2 — 3), 
,3,3 2(n 2 + 3) V ; ' 

c + n 3 

n 3l3 , 3 = 

Integrality of ^1,2,3 requires c = 3n (mod n 2 + 3). If we write c = 3n + a(n 2 + 3) for 
integer a, then we compute <5 2 = (n 2 +3) 2 (a 2 (n 2 +4) + 2an — 3). Hence 5 = (n 2 + 3)/? 
where (3 is integral and (3 2 = a 2 (n 2 + 4) + 2an — 3. Note that in particular this 
forces a 7^ 0. Rewriting it as f3 2 = (an + l) 2 + 4a 2 — 4, we see that (3 = an + 1 
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'mod 2). Thus we compute 

£2 



"1,1,2 

"1,1,3 

"1,2,2 
"1,2,3 
"1,3,3 

™2,2,2 
"2,2,3 



((n 2 + 1)13 - 3n 2 - an(n 2 + 3) - 1), 
ei£2, 



-n(3 + 3n + a(n 2 + 2)), 
(/3 + ara + l), 



2 

1 

£2 

2 

—e±a, 

fw- 

ei£2 



(-/3 + ara + l), 



an 



1), 



((n 2 + + 3n 2 + an(n 2 + 3) + 1), 



^2,3,3 
™3,3,3 



1 

2 

ei£ 2 



n/3 + 3n + a(n 2 + 2)), 
— /3 — an — 1), 



2 

a + n, 

and these are all integral. Nonnegativity of these entries gives further restrictions 
on the parameters. Looking at 71^2,2 + "1,3,3 — ^2/3 we see that e 2 is the sign of (3 
(or (3 = 0, but this gives a = ±1, n = —a, c = a and g(x) = (x — a)(x 2 — 1) which 
is reducible). Looking at 711,1,2 + "2,3,3 = — we see that e\ = —1. Looking 

at 711,2,3 we see that a > 0. Nonnegativity provides additional constraints on the 
parameters, but instead we look at the twist inequalities. 

We saw above that a > hence c = 3n + a(n 2 + 3) > n 2 + 3n + 3 > 0. Thus two 
of the roots d 3 , A 3 ,i and A 3 ,2 of g must be positive. By symmetry, we may assume 
d 3 > ^3,1 > > A 3) 2- Then e 2 = 1 and we have 



S = 



( 1 d 3 /A 3 ,i -d 3 /X 3> 2 d 3 \ 

d 3 /\ 3 ,i 4/A 3 ,2 -1 d 3 

-d 3 /\ 3 ,2 -1 d 3 /\ 3t i -d 3 

V d 3 d 3 -d 3 nd 3 J 



Let M = max(l/A 3 ,i, 1/|A 3)2 |) so that Md 3 = max(d 1 ,d 2 ). Since D 2 = (n 2 + 3)d\ 
and 

- + T^ + — = " 2 + 2, 

A 3,l 



A 3,2 



the diagonal twist inequality coming from the (0, 0) entry gives 

2M 2 <n 2 + 3 + ^4±5. 

d 3 

This inequality allows only finitely many choices of the parameters. 
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If n > 0, then g(—l/<f> n ) = — 0~ 3 < and therefore A 3j 2 > —l/<f> n . Thus M > <p r . 
and we get 

n ,9 9 + 3 

Further since <i 3 A 3j i = Mc and d 3 > A 3; i, we have 



d 3 > \/^c > \J(f) n (n 2 + 3n + 3). 



For n > 2, we get a contradiction by noting that <p n > n, hence these equations 
force 

2n 2 < n 2 + 3 + -=, 

a contradiction. For n — 1, plugging in gives a contradiction. 

If n = 0, then /3 2 = 4a 2 — 3, hence a = 1 and g(x) = x 3 — 3a; 2 + 3. Since 
§3,3 = 0, the (3,3) entry of the twist equation gives 9 2 d\{\ + Q\ + 9 2 ) = 0, hence 
9 1 = 9 2 = e ±2 ™/ 3 . The (0,3) and (1,3) entries give 

3 (1 + - d 2 9 2 ) =D + = 9 1 9 3 (d 1 - d 2 9 ± + 9 2 ). 

This case is realized by (Al, 7)i. 

If n < 0, then 1/M = A 3i i. One easily checks that M and rf 3 are increasing 
functions of c, therefore it suffices to check that the inequality 2M 2 < n 2 + 3 + 
(n 2 + 3) 1//2 <i 3 1 fails for a = 1 and hence c = n 2 + 3n + 3. The inequality fails for 
n < —2. (To see this simply compute both sides for n = —2. For n < —3, note 
that g(-l/n) = -((n + l)/nf < and g((n + l) 2 ) = -(n + l) 4 + n + 2 < 0. 
Therefore M < -n and d 3 > (n + l) 2 . Hence we have 2M 2 > 2n 2 > n 2 + 4 and 
(n 2 + S) 1 / 2 /^ < 1, but these combine to contradict the inequality.) For n = — 1, 
a — 1, the inequality holds, but the resulting polynomial g(x) = x 3 — x 2 — x + 1 
is reducible. Moving up to the next case n — — 1, a — 3 the inequality fails. Thus 
there are no solutions in this case. 

With the cases above completed, we consider G which is not transitive and 
contains no 3-cycle. Up to symmetry, it follows that G must be a subgroup of 
Z/2Z x Z/2Z= ((0 1),(2 3)). 

Case 4. G contains the transposition a = (2 3). 

In this case the parity condition gives eie 2 e 3 = —1. Three instances of the 
usual identity give d 2 = s , 2 = e 2 s 0:3 = e 2 d 3 , d 3 = s , 3 = £3^0,2 = £3^2, and 
d\ = 5o,i = e iSo,i = £idi- Since the di are positive we conclude e\ — e 2 — e 3 = 1, a 
contradiction. Thus we are left with only three possibilities. Either G = Z/2Z = 
((0 1)(2 3)), G = Z/2Z = ((0 1)), or G is trivial. 



Case 5. G contains a — (0 1)(2 3). 
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Using the identities §j : k = ^a(j)^kS a (j),cr- 1 (k) and e\t 2 e 3 = 1 gives 

/ 1 d x d 2 d 3 \ 

g _ d\ 6\ t 2 d 3 e l e 2<^2 
d 2 € 2 d 3 §2,2 §2,3 

\d 3 e x e 2 d 2 §2,3 eiS2,2/ 

Suppose first that e± — 1. Then, orthogonality of the first two columns of 5 forces 
e 2 = — 1 and rfi = d 2 d 3 . Orthogonality of the last column with the first three gives 
the equations 

s 2 , 2 d 2 + s 2i3 d 3 = d 2 (d\ - 1), s 2:2 d 3 + s 2y3 d 2 = d 3 (d\ _ 1)> an d §2,252,3 = -d 2 d 3 . 
Looking at the cases d 2 = d 3 and d 2 7^ d 3 separately, these solve to give 

§2,2 = -1 and s 2 ,3 = d l = d 2 d 3 . 

I 1 d 2 d 3 d 2 d 3 \ 

d 2 d 3 1 -d 3 -d 2 

c?2 ~d 3 -1 rf 2 C?3 

\ rf 3 -d 2 d 2 d 3 -1 J 

This is exactly the S matrix of Case 2 above. Exactly as in that case, we get 
d 2 = 4>m, d 3 = <p n , and di = d 2 d 3 . The iVj matrices and the twist inequalities are 
the same, hence we conclude (m,n) = (0,1), (0,2), or (1,1). (Here we exclude 
m = n = since it gives G trivial.) In the first case, (to, n) = (0, 1), the possible 
twist matrices are given by 9 3 = e ±47 ™/ 5 , 9 2 = ±i, and 9i = 9 2 9 3 . In the second 
case, (m,n) = (0,2), no twist matrix exists. (To see this, note that the (1, 1) and 
(3, 3) entries in the twist equation give 

e 2 1 (<j>l + e 1 + 9 2 (i>l + e 3 ) = D + 

-9 2 M + 9 l + 9 2 <j ) j + 9 3 )=D + . 
Thus 9 X = ±i9 3 . The (0, 1) and (0,3) entries give 

9 1 (l + 9 1 - 9 2 - 9 3 ) = D + 
9 3 (l- 9 1 + 9 2 - 9 3 ) = D + . 

Subtracting these and using the result above gives 9 2 = ±i. Plugging these equa- 
tions into 1 + 9i4>l + 82 + 9 3 4> 2 2 = D+, gives D + = (1 ± i)(l + 9 3 (f>l). Equating 
squared norms gives 9 3 + 9 3 = 1 — <\)\. However 1 — <\>\ = — 2 — 2\pl < —2, so this 
is impossible.) In the third case, (to, n) = (1,1), the possible twist matrices are 
given by 9 2 = e ±4 ™/ 5 , 9 3 = e ±4ni/5 , and 9 1 = 9 2 9 3 . 



Hence 

S = 
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Next consider the case e± = — 1 so 

(I <h 

S = 



di -1 
d 2 e 2 d 3 
\d 3 -e 2 d 2 



d 2 

e 2 d 3 

S 2 ,2 
§2,3 



d 3 ^ 

—e 2 d 2 

§2,3 



By symmetry under interchanging N 2 and A^ 3 , we may assume e 2 = +1. Since a 
is the only nontrivial element of the Galois group, we conclude that 



di 



n, d 2 + — = r, 
d\ 



and d 3 p 



are integers. Hence di = n - If < 0, then d\ < 1 and since it is the largest root 
d 3 < d 2 - But then it follows that trace(iVi) = d 1 — \ jd\ + d 3 /d 2 — d 2 /d 3 < 0, an 
impossibility. Further if n = 0, then the same argument shows d\ — 1 and d 3 = d 2 . 
However, looking at the eigenvalues of N 2 shows cr(d 2 ) = d 3 jd\ > and looking at 
N 3 shows cr(d 3 ) = —d 2 jd\ < 0. Thus n > 1 and d\ is irrational. Since G = Z/2Z, 
it follows that K = Q[di]. Hence we can write d 2 = ad 1 + b and d 3 = dd\ + b for 
rational a, a, 6, 6. Hence 

r = d 2 + -y- = (a + b)d\ + b + a — nb, and s = d 3 p = (a — 6)di + 6 — a + nb. 

d\ d\ 

Hence b = —a, a = b, r = na + 2b and s = nb — 2a. Note that in complex terms 
this gives d 2 + id 3 = (a + ib)(d 1 — i) and r + is = (a + ib)(n — 2i). In particular 
D 2 = 1 + d\ + d\ + d\ = (1 + d?)(n 2 + r 2 + s 2 + 4)/(n 2 + 4). 

§23 = 1 + d\. Since s 2j2 /d 2 
s 2, 2$ 2,3 /{d 2 d 3 ) is an integer. Further 
and hence s 23 = 0. Thus 



Since the columns of S are of equal length s\ 2 



is an eigenvalue of A^, (s 2j2 /d 2 )a(s 2j2 /d 2 
§2,2 7^ 0, since s 2j 2 = would force a{s 2>2 

l + d 2 



> 



> 1. 



2d 2 d 3 d 2 d 3 
The twist inequality coming from the (0, 1) entry of S gives 

(1 + d{yi 2 {\ + a 2 + 6 2 ) 1 / 2 = D <2 + 2- 



.d 2 d 3 



<2 + 



< 2 + 



d\ 
2d 2 d 3 
1 + dl 
l + d 2 
d 1 



l + d 2 
d\ 

(1 + di) 2 



d l 



Rewriting this gives using 



+ s 2 < (n 2 + 4) 



Ad\ + hd\ + Adi + 1 
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j2 ^ n 2 



The twist inequality coming from the (0, 0) entry of S is 2d\ < D + D, hence 

< + 8d 2 + 5di + 1. 

It follows that d\ < 7, hence 1 < n < 6. Together with the bound on r 2 + s 2 
above, this leaves only finitely many possibilities (110 of them, insisting that <i 2 
and d 3 be positive, but finite). 
Since G = Z/2Z, the quantities 

d 2 d 3 r 2 + nrs — s 2 
~dT ~ n 2 + 4 ' 

d 3 d 2 (n 2 + A)rs 

— n + 



d 2 d 3 r 2 + nrs — s 2 

$2,2 §2,3 _ in 2 - 4)r 3 - Y2nr 2 s - 3(n 2 - 4)rs 2 + 4ns 3 
d 2 c?3 (r 2 + s 2 )(r 2 + nrs — s 2 ) 

§2,3 §2,2 4nr 3 + 3(n 2 — 4)r 2 s — 12nrs 2 — (n 2 — 4)s 3 
d 2 ds (r 2 + s 2 )(r 2 + nrs — s 2 ) 

s 2 , 2 s 2 , 3 _ (4 - 3n 2 )(r 4 - 6r 2 s 2 + s 4 ) - 2n(n 2 - 12)rs(r 2 - s 2 ) 
d 2 d 3 (r 2 + s 2 ) 2 (r 2 + nrs — s 2 ) 

are all integers. Only 6 of the 110 examples pass these integrality conditions. 
These are (n,r,s) = (1,2,1), (1,3,-1), (2,2,2), (3,2,3), (4,2,4), and (4,3,1). 
The cases with n = s and r = 2 can be ignored since they give a = and 6=1, 
hence d 2 — 1, d 3 — d±, and s 2j 2 = — 1- Thus s 3i3 = 1. Invoking the symmetry 
under interchanging Ni and puts us back in the case e\ — \. These are just the 
(0, n) examples discussed above. The remaining two examples fail to give integral 
Ni matrices and also fail the twist inequalities. Thus there are no new examples 
in this case. 

Case 6. G contains the transposition a = (0 1). 

Since we can exclude the cases above, a must be the only non-trivial element of 
G. Up to symmetry there are two cases. The parity condition gives eie 2 e 3 = —1. 
Since d 2 = s 0j2 = eie 2 s"i j2 and s 12 = e 2 s 2 = e 2 d 2 , we conclude e\ — 1 and e 3 = — e 2 . 
Then s 2j3 = e 2 e 3 s 3j2 and the fact that S is symmetric forces s 2j3 = s 3j2 = 0. By 
symmetry we may assume e 2 = 1 and e 3 = — 1. Thus we get 

(I di d 2 d 3 \ 
d\ 1 d 2 —ds 

d 2 d 2 §2,2 

\d 3 -d 3 s 3j3 y 
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Note that A2,2 = £2,2/ ^2 and A3 t 3 = 53,3/^3 are integers. Note that orthogonality 
of the third column of S with the first forces A2,2 < 0. Since d\ > a(di) = l/d\, 
we conclude d\ > 1. Thus orthogonality of the first and fourth columns of S gives 
^3,3 > 0. It is straightforward, though somewhat tedious, to build the Ni matrices 
in this case and worry about their integrality and nonnegativity; however, there 
is an easier approach. Since §2,3 = s 3} 2 = 0, the twist equation for the (2, 3) entry 
becomes 

e 2 03{d 2 d 3 - m 2 4) = 0. 

Thus we conclude Q\ = 1. Using this fact the (2,2) entry becomes 

6 2 2 (2di + 62§l 2 ) =D + § 2 ,2- 

Since §2,2/^2 = ^2,2 < is integral and \D + \ 2 = D 2 = (A22 + 2)rf|, equating the 
squared norms of the sides of this equation gives 

02 + 02 = 1-4-- 
*2.2 

The left hand side is an algebraic integer, hence we conclude A2,2 — — 1 and 
62 = exp(dz27ri/3). Similarly, the (3, 3) entry gives A3 5 3 = 1 and 9 3 = exp(±27ri/3). 
Equating the squared lengths of the last two columns of S now gives d 2 = d\ 1 hence 
g?2 = ^3- This is a contradiction, since (7(^2) = A2,i = d 2 jd\ but oid^) = X^i = 
-d 3 /d t . 

Case 7. G is trivial. 

This case is also contained in [CZJ. 

In this case all the di and \j are integral. By symmetry, we may assume 
1 < d\ < d 2 < d 3 . Since every column of S must have squared length D 2 , we see 
that D 2 must be a multiple of d 2 for all i. If d\ = c?2 — ^3, then d\ must divide 
D 2 = 3d 2 + 1. Hence d\ — d% — d 3 — 1. Up to symmetry orthogonality of the 
columns of S forces 



S = 



fl 


1 


1 


M 


1 


-1 


-1 


1 


1 


-1 


1 


-1 


V 


1 


-1 


-v 


fl 


1 


1 


M 


1 


1 


-1 


-1 


1 


-1 


1 


-1 


\1 


-1 


-1 


1 / 



32 



ERIC ROWELL, RICHARD STONG, AND ZHENGHAN WANG 



For the first S matrix, this gives integral iVj matrices 



iVi 



/o 


1 





^ 




(0 





1 




1 











, N 2 = 











1 











1 


1 











\o 





1 






v> 


1 





0/ 



and N-i 











1\ 








1 








1 








Vi 











The possibilities for the corresponding twist matrix are 



T 



or their complex conjugates. (Note that this is the case m 



(1 












(I 








°\ 





i 













i 














1 













-1 





\o 








-i) 




Vo 








V 



n 



of the form 



found in Cases 2 and 5.) For the second 5* matrix, the compatible T matrices are 
listed in Table [2j 

If the di are not all equal, then d\ < d 3 , hence 1 + d\ < d\ and d\ < D 2 < 3d 2 ,. 
Thus we must have D 2 = 2d 2 , i.e., two of the §^3 are zero and the other two are 
±d 3 . Of course s 0)3 = d 3 . Suppose Sj j3 = ±d 3 is the other nonzero entry in the last 
column. Orthogonality of the first and last columns gives d 3 (l ± dj) = 0, hence 
the lower sign is correct and dj = 1. Thus using symmetry we may assume j = 1. 
Orthogonality of the remaining columns with the last column gives s^i = d\ = 1 
and §2,1 = d 2 . Orthogonality of the third column and the first two gives §2,2 = —2 
and 



S 



(\ 
1 

d 2 
\d 3 



1 
1 

d 2 
-d?, 



d-2 
d2 
-2 




-d 3 





However equality of the squared lengths of the last two columns now gives D 2 

d 2 + 2, an impossibility. 



2d 2 



2d 2 + 4 or d 2 3 



□ 



5. Realization of fusion rules and Classification of MTCs 

In this section, for each modular fusion rule (X; S) in Theorems 13. 1\ \3.2\ 14.14 
we will first determine all modular symbols with this fusion rule which also satisfy 
(1)(2) of Proposition 12. 3( then classify all MTCs realizing each such modular 
symbol. 

For each modular fusion rule (N; S) , there are two choices of compatible S 
matrices: S = j^S or —j^S. When the two modular symbols are realized by 
(2+l)-TQFTs, respectively, one TQFT is obtained from the other by tensoring 
the trivial theory with S = (—1). The quantum invariant of the 3-sphere will be 



ON CLASSIFICATION OF MODULAR CATEGORIES 



33 



jy or —j), respectively. Also the topological central charge c of the two theories 
will differ by 4. Another symmetry for modular symbols is complex conjugation: 
to change (N; S, T) to (N; S\ T'). Complex conjugation of a modular symbol gives 
rise to a different modular symbol if one of the S, T is not a real matrix. 

Given an S matrix, we can obtain all fusion matrices by using the Verlinde 
formulas. Instead of listing fusion matrices, we will present them as fusion rules. 
In the following, we will not list trivial fusion rules such as l<S>x = x and those that 
can be obtained from obvious identities such as x ®y — y ® x. We will also write 
x<g)y as xy sometimes. Then we use relations (3)(i)-(iv) of Definition 12. II together 
with (1) (2) of Proposition 12.31 to determine the possible T-matrices. As Example 
12.21 illustrates, Proposition 12.31 is necessary to get finitely many solutions in some 
cases. We find that there are finitely many modular symbols (N; S, T) of rank< 4 
satisfying Proposition ^. 31 Modulo the symmetry S — > —S, these modular symbols 
are classified in Table [21 In the table, C, m = e~^\ The labels will be {1,X} for 
rank=2, {1, X, Y} for rank=3, and {1, X, Y, Z} for rank=4. They will correspond 
to rows 1,2,3,4 of the S matrices. The # is the number of modular symbols 
satisfying Proposition 12.31 modulo the symmetry S — ► —S. The column P stands 
for primality, and the column G is the Galois group of the modular fusion rule. 

With modular symbols determined, we turn to realizing each of them with 
MTCs. First let us consider the S-matrices corresponding to Theorem !4.1( 3)-(5). 
In each of these cases there is a rank=2 tensor subcategory corresponding to the 
objects labelling columns 1 and 3 of the S matrix. Further inspection shows that 
the submatrix of S corresponding to rows and columns 1 and 3 is invertible. It 
is obvious that the tensor subcategory generated by the trivial object and the ob- 
ject labelling column 3 is a modular subcategory. For the S'-matrices of Theorem 
14. 1( 3). (4) these rank=2 modular subcategories are equivalent to the UMTCs cor- 
responding to Theorem 13.1( 1). while the modular subcategory corresponding to 
the S'-matrix of Theorem 14.1( 5) is equivalent to (one of) the UMTCs coming from 
Theorem 13.1( 2). By [M2j [Theorem 4.2] this implies that the MTCs correspond- 
ing to these S matrices are direct products of rank=2 MTCs. For this reason we 
will not write down realizations or complete data for these MTCs as they can be 
deduced from their product structure. 

MTCs realizing the remaining 8 nontrivial modular symbols are prime, i.e. they 
do not have non-trivial modular subcategories. To complete the classification, 
we need to solve the pentagon and hexagon equations for all 8 modular symbols. 
The solutions of the pentagon equations are organized into the F-matrices whose 
entries are called 6j symbols. The solutions of hexagons are given by the braiding 
eigenvalues. 

5.1. F-matrices. Given an MTC C. A 4-punctured sphere S^ bcd , where the 4 
punctures are labelled by a, b, c, d, can be divided into two pairs of pants(=3- 
punctured spheres) in two different ways. In the following figure, the 4-punctured 
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Table 2. Rank< 4 unitary modular symbols 



S matrix 


Fusion rules 


T matrix 


#s 


P 


G 


S = (l) 




(1) 


1 


Yes 


1 


Thm 3.1(1) 


X 2 = 1 


Diag(l,±i) 


2 


Yes 


1 


Thm 3.1(2) 


X 1 = 1 + X 


Diag(l, (( 5 ) ±2 ) 


2 


Yes 


Z 2 


Thm 3.2(1) 


X 1 = X*, XX* = 1, 
(X*) 2 = X 


Diag(l,C 3 ±1 ,C3 ±1 ) 


2 


Yes 


Z 2 


Thm 3.2(2) 


X 2 = 1 + Y, 
XY = X, Y 2 = 1 


Diagfl, (Cie) 2fc+1 , -1) 


8 


Yes 


z 2 


Thm 3.2(3) 


X 2 = 1 + X + Y, 
XY = X + Y, 
Y 2 = 1 + X 


Diag(l, (Cr) ±5 , (C 7 ) ±x ) 


2 


Yes 


z 3 


Thm 4.1(1) 


X 2 = Y = (X*) 2 , 
XX* = 1 = Y 2 , 
XY = X*, X*Y = X 


Diag(l,-l,(C 8 ) ±m , (Cs) ±m ), 
m=l,3 


4 


Yes 


z 2 


Thm 4.1(2) 


X 2 = 1, XY = Z, 
XZ = Y,Y 2 = 1, 
YZ = X, Z 2 = 1 


Diag(l, -1,6!, ex), 
e? = l 


2 


Yes 


1 


Thm 4.1(3) 


X 2 = 1, XY = Z, 
XZ = Y,Y 2 = 1, 
YZ = X, Z 2 = 1 


Diag(l, ^,02,^2), 
= -1 


3 


No 


1 


Thm 4.1(4) 


X 2 = 1 + X, 
XY = Z, 

XZ — Y + Z,Y 2 — 1, 
YZ = X, Z 2 = 1 + X 


Diag(l, ^^W, 

01 = (C 5 ) ±2 , 02 = ±1 


4 


No 


z 2 


Thm 4.1(5) 


X 2 = 1 + X, 

xr = z, 

XZ = F + Z, 

r 2 = i + y, 

YZ = X + Z, 

Z 2 = l + X + Y + Z 


Dia g (i, e 1 ,e 2 ,e 1 e 2 ), 
0i = (C 5 ) ±2 , 02 = (C 5 ) ±2 


3 


No 


z 2 


Thm 4.1(6) 


x 2 = i + x + y, 
xy = x + y + z, 

XZ = Y + Z, 

Y 2 = l+X + Y + Z, 

yz = x + y, 

Z 2 = 1 + X 


Diag(l,C 9 ±2 ,(C 9 ) ±6 ,(C9) T6 ) 


2 


Yes 


z 3 
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sphere is the boundary of a thickened neighborhood of the graph in either side, and 
the two graphs encode the two different pants-decompositions of the 4-punctured 
sphere. The F-move is just the change of the two pants-decompositions. 

When bases of all pair of pants spaces Hom(a <g> b, c) are chosen, then the two 
pants decompositions of S^ bcd determine bases of the vector spaces Hom((a <g> b) ® 
c, d), and Hom(aCg) (b<S>c), d), respectively. Therefore the F-move induces a matrix 
F%' b,c : Hom((a®6) <S>c,d) — > Hom(a® (&<g)c), d), which are called the F- matrices. 
Consistency of the F matrices are given by the pentagon equations. 

For each quadruple (a,b,c,d), we have an F-matrix whose entries are indexed 
by a pair of triples ((m,s,t), (n,u,v)), where m,n are the labels for the internal 
edges, and s, t, u, v are indices for a basis of the Kom(x®y, z) spaces with dim > 1. 
For the MTCs in our paper, none of the Hom(a ® b, c) has dim > 1, so we will 
drop the s, t, u, v from our notation. If one of the a, b, c in F^' b,c is the trivial label, 
then we may assume F%' ,c is the identity matrix. But we cannot always do so if 
d is the trivial label. In the following, any unlisted F matrix is the identity. 



a b c a b c 




d d 



5.2. Braidings and twists. The twist of the simple type X { will be denoted by 
$i, and it is defined by the following positive twist. 



The braiding eigenvalues are defined by the following diagram: 



a b a b 




c c 

The consistency equations of the braidings are given by two independent families 
of hexagon equations. If c is the trivial label, and label a is self-dual, then R® a = 
v a Q~a^ , where v a is the Frobenius-Schur indicator of a. 
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5.3. Explicit data. In this section, we give the explicit data for at least one 
realization of each prime modular fusion rule. Since each modular symbol can 
have up to 4 MTC realizations, we will present the complete data for only one of 
them. We choose one with the following properties: 

(1) The (0, 0) entry of the S matrix is j=j, where D is the total quantum order; 

(2) In a category with a generating non-abelian simple object X, we choose 
a theory with the positive exponent of the twist 9x being the smallest. This is 
inspired by anyon theory that the simple object with the smallest exponent is the 
most relevant in physical experiments. 

If a modular symbol (N; S, T) is realized by an MTC, then the modular symbol 
(N;—S,T) is also realized by an MTC. The modular symbol (N; S', T') is real- 
ized by complex conjugating all F matrices and braidings of (N; S,T). So in the 
following each group of data will be for 4 MTCs if any of S, T, F and braidings 
are not real; otherwise there will be two. We choose the F matrices to be unitary, 
and real if possible. 

In anyon theory, labels will be called anyon types. The smallest positive expo- 
nent of a twist 0i will be called the topological spin of the anyon type i. Topological 
spins are the conformal dimensions modulo integers of the corresponding primary 
field if the MTC has an corresponding RCFT. The last line of the data lists all 
quantum group realizations of the same theory. We did not list the Frobenius-Schur 
indicators of anyons because they can be calculated by the formula in Proposition 
12. 31 In the following data, only the semion s and the (Ai, 2) non-abelian anyon a 
have Frobenius-Schur indicator=— 1. 

5.3.1. Semion MTC. We will use s to denote the non-trivial label. 
Anyon types: 

Fusion rules: s 2 = 1 
Quantum dimensions: {1,1} 
Twists: 8i = 1, 9 S = i 
Total quantum order: D = \[2 
Topological central charge: c = 1 
Braidings: Rf = i 

S-matrix: S = -j= 

F-matrices: F!!' s ' s = (— 1) 
Realizations: (A±, 1), (Ej, 1). 

5.3.2. Fibonacci MTC. We will use (f to denote the golden ratio ip 
the non-trivial label. 

Anyon types: {1,t} 
Fusion rules: r 2 = 1 + r 
Quantum dimensions: {1,^} 




= ^fl and r 
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4-7TZ 



Twists: e 1 = l,9 T = e— 
Total quantum order: D = 2cos(^) = 2s h^^) 
Topological central charge: c = ^ 
Braidings: K[ — e 5 , it/ = e 5 

S-matrix: 5 = ^ ^ ^ , 
F-matrices: F T ™" = ^) 

Realizations: (Ai,3)i, (G 2 , 1), complex conjugate of (F 4 , 1). 

5.3.3. Z 3 MTC. We will use 00 for both a non-trivial label and the root of unity 
to = e 27 ™/ 3 . No confusions should arise. 
Anyon types: {l,cu,ui*} 



Fusion rules: u 2 = u*, uu* = 1, (u*) 2 
Quantum dimensions: {1, 1, 1} 
Twists: Q\ — 1, 9u — Ou* — 
Total quantum order: D = \/3 
Topological central charge: c = 2 

Braidings: R^ = Rf>" = e"¥ ?J R^ = Eg^ = e~^, 

'i 1 r 

S-matrix: 5 = 4^ | 1 c<j c<j 2 



1 cj 2 UJ 
- (1) for a 
Realizations: (A 2 ,l), (E 6 ,l) 



F-matrices: F^' b ' c = (1) for any a, b, c, d, 



5.3.4. Ising MTC. We will use l,cr,ip to denote the non-trivial labels. 
Anyon types: {1,(7, ip} 

Fusion rules: a 2 = 1 + if), aif) = ipa = a, if) 2 = 1 
Quantum dimensions: {1, y/2, 1} 
Twists: #i = 1, 9 a — eir, 0^, = — 1 
Total quantum order: D = 2 
Topological central charge: c = | 

Braidings: i2f ff = e~f , = -1, i?^ 7 = = -i, R% = e¥ 
( 1 V2 1 

-V2 I, 

V 1 -V2 1 



S-matrix: 5 = \ 



F-matrices: = ^ (\ _\) = (-l),F^' a = (-1), 

Realizations: complex conjugate of (E 8 ,2). 
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5.3.5. (Ax, 2) MTC. We will use l,a, if) to denote the non-trivial labels again. 
Anyon types: {l,a,ip} 

Fusion rules: a 2 = 1 + ip, = V' ' = cr,ip 2 = 1 
Quantum dimensions: {1, ^/2, 1} 
Twists: Q\ = 1, 9 a = e~8~, 6^ = — 1 
Total quantum order: D = 2 
Topological central charge: c = | 



Braidings: Rf = -e"¥, Rf* = -1, flf" = Rf = i, R™ 
( 1 
V2 



S-matrix: S = \ 



Realizations: (Ax, 2). 



V2 









-V2 








V2 [l 





5.3.6. (Ax,5)i MTC. We will use l,a, (3 to denote the non-trivial labels. Note 
that l,a,/3 are special labels for 1,Y,X in Theorem 3.2(3) of Table [2j 
Anyon types: {l,a,j3} 

Fusion rules: a 2 = 1 + (3, a(3 = a + (3 1 [3 2 = l + a + (3 
Quantum dimensions: {1, rf, d 2 — 1}, where d = 2 cos(y) 

Twists: 8x = 1, 9 a = e¥ , 9p = 
Total quantum order: D = 2si ^J^ 
Topological central charge: c = y 
Braidings: Rf = e -2 ^,/?^ = e _i ¥, 

Rf = e~^,R a / = R? a = e¥, R*> = e ~¥ 

/ 1 d rf 2 - l\ 

S-matrix: S = -rf 2 + 1 1 

\rf 2 -l 1 -rf / 

F-matrices: see the end of this subsection. 
Realizations: (Al,5)i 



5.3.7. Z4 MTC. We will use l,e, cr, a* to denote the non-trivial labels. 
Anyon types: {1, e, a, a*} 

Fusion rule: e 2 = 00* = 1, a 2 = (a*) 2 = e,ae = a* , a*e = a 
Quantum dimensions: {1, 1, 1, 1} 
Twists: 61 — 1, 9 e — —1, 9 a = 9 a * = 
Total quantum order: D = 2 
Topological central charge: c = 1 
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(7 ,{J* 



Braidings: TCf = -1,R^ = R 

R a*,e = R e,o* = _ { 

(11 1 l\ 



2LL T)G,G* 

e 4 , R{ 



e 4 , R a ', 



S-matrix: S — \ 



11-1-1 
1 -1 -i i 
\1 -1 % -ij 
F-matrices: F°f> c = Ff> a *> a * = F^ ( - IV 
Realizations: (A 3 , 1), (D 9 , 1). 



(T* 



F a*,e,a = (_1) 



5.3.8. Tone code MTC. The fusion rules are the same as Z 2 x Z 2 , but the theory 
is not a direct product. We will use l,e,m,e to denote the non-trivial labels. 
Anyon types: {1, e, m, e} 

Fusion rules: e 2 = m 2 = e 2 = 1, em = e, ee = m, me = e 
Quantum dimensions: {1, 1, 1, 1} 

Twists: e 1 = e e = e m = i, e e = -1 

Total quantum order: D = 2 
Topological central charge: c = 
Braidings: i?i' e = -l,Rf m = l,R t 

~ l i R ni = l i R m = _1 

/ll 1 1 \ 
11-1-1 



m,e 
e 



1 D e > e — z? m 
1, rt^ — iti 



l,R; 



e,m 
e 



1, RT' e = 



S-matrix: S — h 



1 
V 



-1 1 -1 
-1-1 1/ 

F-matrices: F%' b ' c = (1) for all a, 6, c, d. 
Realizations: (D 8 , 1), D(Z 2 ) — quantum double of Z 2 . 



5.3.9. (_D 4 , 1) MTC. The fusion rules are the same as Z 2 x Z 2 , but the theory is 
not a direct product. We will use l,e,m, e to denote the non-trivial labels again. 
Anyon types: {1, e, m, e} 



Fusion rules: e 2 



m 



1, em = e, ee = m, me = e 



Quantum dimensions: {1, 1, 1, 1} 

Twists: e l = i,e e = e m = e e = -l 

Total quantum order: D = 2 
Topological central charge: c = 4 
Braidings: R{' e = -l,R e e ' m = 
1, R™' e = -1, = 1, R^ = -1 

/ll 1 1 \ 
11-1-1 



m,e 



1, iff' 



■1 R e ' m = 

L i - LL e 



S-matrix: S — k 



F-matrices: F d ' ' 
Realizations: (_D 4 ,1). 



1-11-1 
VI -1-1 I J 
- (1) for all a, b, c, d. 
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5.3.10. (Ai, 7) l MTC. We will use 1, a, u, p to denote the non-trivial labels. Note 
that 1, a, u, p are special labels for 1, Z, Y, X in Theorem 4.1(6) of Table [2] 
Anyon types: {l,a,u,p} 

Fusion rules: a 2 = 1 + u, au = a + p, ap = oj + p, J 2 = 1 + uj + p, 
up = a + U + p, p 2 = l + a + uj + p 

Quantum dimensions: {1, d, d 2 — 1, d + 1}, where d — 2 cos(|) and d 3 = 3d + 1. 
lwists: v\ — 1, v a — e a , W = e » ,v p = e s 
Total quantum order: D = . 3 , 7r , 

I 2sm(|) 



Topological central charge: c = y 
Braidings: i?? Q = e" 2 ^, iif = e - ^, i*f = 
C = i?r = e™, = ^ = iC = -1 

i? 2 P = iC = e^, iC = = e-¥, = e¥, Rpp = e~f,R^ = e 2 ? , 
R au = R ua = e -^ jR ap = R P» = e ~f,R^ = = e*? , Rf = e^,R p p p 

2tvi 



S-matrix: S = 



(I d d 2 - 1 d+1 \ 

d -d-1 d 2 -l -1 
d 2 - 1 d 2 - 1 -d 2 + 1 
-1 -d 2 + l d j 
F-matrices: see below. 

Realizations: (Ax, 7)i, complex conjugate of (6*2)2- 

The list of all F matrices for an MTC can occupy many pages. But they are 
needed for the computation of quantum invariants using graph recouplings, the 
Hamiltonian formulation of MTCs as in [LWe] or the study of anyon chains [FTLj . 
For the MTCs (Ai,k)i with odd k, all the data of the theory can be obtained 

from |KLj . For k = 5, choose A = ie~^& , the label set is £ = {0, 2, 4} in [KLj and 
= 1,4 = a, 2 = /3. For k = 7, set A = ieW, the label set is £ = {0,2,4,6} in 
[KL] and = 1, 6 = a, 2 = w,4 = p. The twist is given by 9 a = (-\) a A a ^ a+2 \ and 

the braiding it™ = (—1) 2 A 2 . lhe formulas tor 6j symbols can 

be found in Chapter 10 of [KLj . The F matrices from [KLJ are not unitary, but the 
complete data can be presented over an abelian Galois extension of Q. To have 



unitary F matrices, we need to normalize the 6 symbols as 9(i,j, k) = ^Jdidjdk- 

The (Ai, k)i,k odd, MTCs have peculiar properties regarding the relation be- 
tween the bulk (2 + l)-TQFTs and the boundary RCFTs. To realize (A±, k)i, k 

odd, using the Kauffman bracket formalism, we set A = ie 4 ( fc + 2 ) . In order to fol- 
low the convention above, we choose A = ie 4 ( fc + 2 ) if k = 1 mod 4, and A = ie 4 ( k + 2 '> 
if k — — 1 mod 4. Note that in both cases A is a 2(k + 2)th root of unity. We have 



(5.1) (Ai, k) = (Ai, k)i x the semion, 
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if k 



1 mod 4, and 



(5.2) 



(Ai, k) = (Ai, k)i x the semion, 



if k = — 1 mod 4. The central charge of (Ai,k) is -M^, which implies that the 
central charge of (Ax, k) i is Ck = 1 — if & = 1 mod 4, and = 1 + if 
= — 1 mod 4. 

In Table El we list all unitary quantum groups categories of rank< 12 from the 
standard construction. For notation, see |HRWj . 

Remark 5.1. The following serves as a guide to Table [3j 

(1) In general we will list these categories as (X r , k) for the category obtained 
from a quantum group of type X r at level k. Observe that the correspond- 
ing root of unity is of order I = mk + h where m — 1 for X = A, D, or E; 
m = 2 for X = B, C or F and m = 3 for X = G, and h is the dual Coxeter 
number. 

(2) The category (A r ,k) has a modular subcategory (A ri k)^_ generated by 

r + l 

the objects with integer weights provided gcd(r + 1,A;) = 1. These are 
found on line 5 of Table [3] where 



L = {(1, 2s + 1), (2s, 2), (2, 4), (2, 5), (2, 7), (3, 3), (4, 3), (6, 3) : 1 < s < 11}. 



(3) We include the examples of pseudo-unitary categories coming from low- 
rank coincidences for quantum groups of types F 4 and G 2 at roots of unity 
of order coprime to 2 and 3 respectively. 

(4) This list includes different realizations of equivalent categories. We elim- 
inate those coincidences that occur because of Lie algebra isomorphisms 
such as sp 4 = SO5 etc., and do not include the trivial rank= 1 category. 

(5) NSD means the category contains non-self-dual objects. 

(6) "c.f. (X r , k)" means the categories in question has the same fusion rules 
as those of (X r , k). 

(7) We include the three categories coming from doubles of finite groups with 
rank<12, although they are not strictly speaking of quantum group type. 

5.4. Classification. In this section, we explain Table [TJ We identify MTCs whose 
label sets differ by permutations. For the trivial MTCs, the two MTCs are distin- 
guished by the S matrices: S — (±1). 

For the Z 2 fusion rule, unitary MTCs are the semion MTC and those from the 
two symmetries S — > — S and complex conjugate. 

For the Fibonacci fusion rule, unitary MTCs are the Fibonacci MTC and those 
from the two symmetries S — > —S and complex conjugate. 

For the Z3 fusion rule, all unitary MTCs are the one listed in last subsection 
and those from the two symmetries S — ► — S and complex conjugate. 



42 ERIC ROWELL, RICHARD STONG, AND ZHENGHAN WANG 

Table 3. Unitary Quantum Group Categories of rank < 12 



(X r , k) 


Rank 


Notes 




(A., 1), r < 11 


r + 1 


r > 2 NSD, abelian 


r + 2 


{M, k), k < 11 


k + 1 




k + 2 




6 


NSD 


5 




10 


NSD 


6 


(A 1c\ i (r h\ C T 

r + 1 


1 (k+r\ 
r+1 V k ) 




ft, -|- / -|- 1 


(-Br, 1) 


3 


c.f. (Al,2) 


4r 


(S r ,2), r < 8 


r + 4 


nnite braid imager 


4r + 2 


(#2,3) 


10 




12 


(C r , 1) r < 11 


r + 1 


c.f. (Ai,r) 


2(r + 2) 


(C 3 , 2) 


10 




12 


(£>2r, 1) 


4 


r even c.f. D UJ (Z 2 ) 


4r — 1 


(D 2r +i, 1) 


4 


C.f. (A3, 1) 


4r + 1 


(A-,2), r = 4,5 


11, 12 


r = 5 NSD 


8, 10 


/ 771 1 \ 1 1 O 

{E 6 ,k), fc = 1,2 


3, 9 


NSD 


13, 14 


(£■7, k), 1 < fc < 3 


2,6,11 




19,20,21 


(£ 8 ,fc), 2<A;<4 


3,5,10 




32,33,34 


(F 4 , fc), 1 < k < 3 


2,5,9 




20,22,24 


(G 2 ,A;), 1< fc < 5 


2,4,6,9,12 




15,18,21,24,27 




10 


c.f. (£ 8 ,4) 


17 


G 2 


5,8,10 




11,13,14 


D"(Z 2 ) 


4 


prime 




D»(Z 3 ) 


9 


prime 




D«(S 3 ) 


8 


c.f. (£4,2) 





For the Ising fusion rule, there are a total of 16 theories divided into two groups 
according to the Frobenius-Schur indicator of the non-abelian anyon X, X 2 = 
1 + Y. There are 8 unitary MTCs with Frobenius-Schur indicator=l. Their twists 
are given by 6x = e~ for m — 1,7, 9, 15. The Ising MTC is the simplest one with 
m = 1 and central charge c—\. The theory m = l,m— 15 are complex conjugate 
of each other, so are the m = 7, 9. The other 4 MTCs are obtained by choosing 
—S. There are 8 unitary MTCs with Frobenius-Schur indicator^— 1. Their twists 
are 9 X = e^f 1 for m — 3,5, 11, 13. The SU(2) at level k = 2 is the simplest one 
with m = 3 and central charge c = |. The MTCs m = 3 and m — 13 are complex 
conjugate, so are m — 5, 11. The other 4 are those with —S. The Ising MTC is 
not an SU{2) theory. It can, however, be obtained as a quantum group category 
as the complex conjugate of E 8 at level=2. Note that the F matrices in each group 
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of 8 are the same, but their braidings are different. The SU(2) level=2 theory has 
Fx XX — —F° aa with the other F matrices the same as the Ising theory. 

For the {A\, 5) 1 fusion rule, all unitary MTCs are the one listed in last subsection 
and those from the two symmetries S — > — S and complex conjugate. 

For the Z2 x Z2 fusion rules, there are two groups of theories depending whether 
or not the theory is a product. There are 4 theories which are not direct products, 
and 6 product theories. The toric code MTC has another version, which could 
also be called the toric code: it has 9 e = 9 m = —1. All F matrices are 1. The 
braidings Rf e = R™ m = Rf™ = — 1, R™ e = 1, and others are the same as the toric 
code. Another two are the —S versions. The product theories are the products 
of the semion MTC and its complex conjugate. There are 4 possible theories, but 
two of them are the same: semion x complex conjugate is the same as complex 
conjugate x semion. Hence there are 3 theories here. With the — S versions, we 
have 6 product theories. 

For the Z 4 fusion rule, the Galois group action of the MTC listed above is Z 4 . 
Its actions give rise to 4 theories with 6x = Ox* = for m — 1,3, 5, 7. They all 
have the same F matrices. The —S versions give a total of 8. 

For the (Ai,3) fusion rule, this is the product of the semion fusion rule with 
the Fibonacci fusion rule. There are 4 product theories from semion, Fibonacci 
and their complex conjugates. These 4 theories are different, and the other 4 
come from their —S versions. Let us choose the product of the semion with the 
Fibonacci as a representative theory, then we have 4 anyons, l,tp,r,s, where r is 
the Fibonacci anyon, and ip is the same as r tensoring the semion s. 

For the (A\, 7) 1 fusion rule, all unitary MTCs are the one listed in last subsection 
and those from the two symmetries S — > — S and complex conjugate. 

The analysis of the Fibonacci x Fibonacci fusion rule is the same as that of the 
semion x semion fusion rule. 

6. Conjectures and Further Results 

In this section we briefly discuss several conjectures concerning the structure 
and application of MTCs. 

6.1. Fusion Rules and the Finiteness Conjecture. Since topological phases 
of matter are discrete in the space of theories, therefore, MTCs, encoding the 
universal properties of topological phases of matter, should also be discrete. 
It is conjectured |Waj : 

Conjecture 6.1. If the rank of MTCs is fixed, then there are only finitely many 
equivalence classes of MTCs. 

By Ocneanu rigidity, this is equivalent to there are only finitely many modular 
fusion rules realizing by MTCs of a fixed rank. 
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Proposition 6.2. There are only finitely many equivalence classes of unitary 
MTCs with total quantum order D < c, where c is any given universal constant. 

Proof. For a unitary rank=n MTC, all quantum dimensions d r > l,r G £. So 
D > y/n. If D < c, then n < c 2 . By Verlinde formula 12.31 we have n* • = 

I S"=o Sir %f kr I < ^ E"=o j; <nD < c 3 for any z,j, fc. Therefore, there are only 
finitely many possible fusion rules. By Ocneanu rigidity, there are only finitely 
many possible MTCs. 

□ 

6.2. Topological Qubit Liquids and the Fault-tolerance Conjecture. Topo- 
logical phases of matter are quantum liquids such as the electron liquids exhibiting 
the FQHE, whose topological properties emerged from microscopic degrees of free- 
dom. This inspires the following discussion. 

Let A be a triangulation of a closed surface E, Ta be its dual triangulation: 
vertices are centers of the triangles in A, and two vertices are connected by an 
edge if and only if the corresponding triangles of A share an edge. The dual 
triangulation of A is a celluation of £ whose 1-skeleton is a tri-valent graph. It 
is well-known that any two triangulations of the same surface S can be transformed 
from one to the other by a finite sequence of two moves and their inverses: the 
subdivision of a triangle into 3 new triangles; and the diagonal flip of two adjacent 
triangles that share an edge (=the diagonal). Dualizing the triangulations into 
celluations, the two moves become the inflation of a vertex to a triangle and the 
F move. 

Definition 6.3. (1) Given an integer k > 0, a A;-local, or just local, qubit 
model on (S, Ta) is a pair (!Ka,-£^a), where JCa is the Hilbert space 
®eer A ^-- 2 5 an d -f^A is a A; local Hamiltonian in the following sense: if a is a 
sum of Hermitian operators of the form id® ■ ■ ■ ®id®> Ok®id® ■ ■ ■ ® id, 
where acts on < k qubits. 
(2) A modular functor V is realized by a topological qubit liquid if there is a 
sequence of triangulations {A i }^ :1 of £ whose meshes — > as i — > oo, an 
integer k, and uniform local qubit models on (E,rAj such that 

(i) the groundstates manifold of each is canonically isomorphic to 
the modular functor V^(E) as Hilbert spaces; 

(ii) the mapping class group acts as unitary transformations compatibly; 

(iii) there is a spectral gap in the following sense: if the eigenvalues of 
the Hamiltonians are normalized such that = Aq < < • • • , then 
X\> c for all i, where c > is some universal constant. 

The scheme for the local qubit models should be independent of the geometry 
of the surface S, and have a uniform local description. The modular functor 
determines a unique topological inner product on V(£). We require that the 
restricted inner products from IHa, to the groundstates of agree with the 
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topological inner product on V(E). To identify the Hilbert space "K^ of one 
triangulation with another, we consider the two basic moves: F move and inflation 
of a vertex. The F move does not change the number of qubits, so the two Hilbert 
spaces !Ka 4 have the same number of qubits. We require that the identification 
to be an isometric. For the inflation of a vertex, the inflated celluation has 3 new 
qubits, so we need to choose a homothetic embedding with a universal homothecy 
constant. 

The action of the mapping class group is defined as follows: consider the moduli 
space of all triangulations of E that two triangulations are equivalent if there dual 
graphs Ta are isomorphic as abstract graphs. By a sequence of diagonal flips, we 
can realize a Dehn twist. Each diagonal flip is an F move, and their composition 
is the unitary transformation associated to the Dehn twist. 

Conjecture 6.4. (1) Every doubled MTC C can be realized as a topological 
qubit liquid. 

(2) The groundstates V(E) = C CKa^ form an error- correction code for 
each triangulation Aj. 

6.3. Topological Quantum Compiling and the Universality Conjecture. 

Every unitary MTC gives rise to anyonic models of quantum computers as in 
|FKLWj . Quantum gates are realized by the braiding matrices of anyons, i.e. the 
afforded representations of the braid groups. Topological quantum compiling is 
the question of realizing desired unitary transformations by braiding matrices in 
quantum algorithms, in particular for those algorithms which are first described 
in the quantum circuit model such as Shor's famous factoring algorithm. 

To choose a computational subspace, we will use the so-called conformal block 
basis for the Hilbert space V(D 2 , a»; aoo) of a punctured disk, where Oqo labels 
the boundary. Conformal block basis is in one-one correspondence to admissible 
labelings m, n, ■ ■ ■ ,p of the internal edges of the following graph subject to the 
fusion rules at each trivalent vertex. As explained in Section I5.1[ the tri-valent 
vertices also need to be indexed if multiplicities n\ - > 1. 

di 02 «3 a m 




The braiding of two anyons Oj, a^+i in a conformal block basis state is represented 
by the stacking the braid on top of the above graph at i, i + 1 positions. 

Definition 6.5. An MTC C has property F if for every object X in C and every 
m the representation of !B m on V(D 2 , X, ■ ■ ■ ,X; a^) factors over a finite group 
for any o m G £. 



46 



ERIC ROWELL, RICHARD STONG, AND ZHENGHAN WANG 

Table 4. Unitary prime MTCs rank<4 



Realization 


PSL(2,Z), Relations 


Property Fl 


Universal Anyons 


XT -i- 

Vectc 


1, o = 1 =1 


Yes 




f A 1 \ 


PbL(2, 3), 1=1 


Yes 






PSL(2,5), T b = I 


No 


T 


(^2, 1) 


PSL(2,3), T 6 = I 


Yes 




(Ai,2) 


PSL(2,8), 

T 16 = ( T 2 5T )3 = j 


Yes 




(^i,5)i 


PSL(2,7), 

T 7 = (T 4 ST 4 S) 2 = I 


No 


a,j3 


(4>, 1) 


PSL(2,8), 

T 8 = (T 2 ST) 3 = I 


Yes 




D(Z 2 ) 


PSL(2,2), T 2 = I 


Yes 




(^i,7)i 


PSL(2,9), 

T 9 = (T 4 ST 5 S) 2 = I 


No 


a, oo, p 



The following is conjectured by the first author (see |NR] ) : 

Conjecture 6.6. Let G be an MTC. 

(a) IfQ is unitary, then it has property F if and only if(di) 2 G N for each simple 
object Xi or, equivalently, if and only if the global quantum dimension 

D 2 en. 

(b) In general, Q has property F if and only if (FPdim(Xj)) 2 G N for each 



simple object Xi, where FPdim is the Frobenius- Perron dimension, i.e. 
Frobenius-P erron eigenvalue of the fusion matrix N^. 



the 



The verification of this conjecture for UMTCs of rank< 4 is summarized in 
Table H 



Theorem 6.7. The following anyons are universal in the sense of [F KLW] : 
the Fibonacci anyon r, the (A\, 5) i anyons a, (3, the {A\, 7) i anyons a, uj, p, the 

two anyons ip,r in (Al,3) (see \5.4\ for notation), and the two t's in Fib x Fib. 
Universal anyonic quantum computation can also be achieved with the anyon 

r xt in Fib x Fib, but images of the representations of the braid groups from this 

anyon are not as large as possible. 

Anyons that correspond to ip,r,a, /3,u), p in other versions are also universal. 



Proof. We deduce the proof from [FLW2j [LRW] [LWa] . 

Universality of tp and r is given in }FLW2j . The anyons a are both the fun- 
damental representations of (Ai,k) up to abelian anyons. The universality of 
fundamental representation anyons are established in |FLW2] . Therefore, both 
oj's are universal. 
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To prove that /3,u, p are universal, we first show that their braid representations 
are irreducible. By inspecting the braiding eigenvalues in Section 15. 3^ we see that 
they satisfy the conditions of |TW] [Lemma 5.5] [HRWJ [Proposition 6.1]. It follows 
that the braid representations are irreducible. Universality now can be proved 
following |FLW2j or [LKW] . 

□ 

Appendix A. Non-self dual rank< 4 MTCs with S. Belinschi 
Every rank=l, 2 MTC is self-dual, so we will start with rank=3. 

A.l. Nonselfdual Rank=3. The three labels will be 0, 1, 2 such that 6 = 0, 1 = 
2, 2 = 1. The modular S matrix is of the form: 




S22 = Sn, §12 = Su follows from sj . = Sjj. Unitarity of S implies 
(A.l) l + d 2 = 2\x\ 2 , 



(A.2) 



d 2 + x 2 + x 2 



0. 



(A.3) 



1 + x + x = 0. 



The fusion matrix Ni has eigenvalues d, §, Their sum d + = d — \ is an 

integer. Their product H- = ^f- = \(\ + d) is also an integer. Therefore, d is an 
integral multiple of |, so d is an integer. 

Let 6 be the twist of label 1, hence of label 2. Using identity ( 12. 21) . we get 

72 , + 0-l 



(A.4) 



2d z 



0. 



Therefore, 2d 2 < 3. Since d ^ 0, the only possible integers are <i 2 = 1, hence 



x 



1. Then 1 + x + x = leads to x = e 3 



A.2. Nonselfdual rank=4. Now we turn to the non-self dual rank=4 case. The 

4 labels will be denoted as 1,Y, X, X* , where Y is self dual and X, X* dual to 

each other. Taking into account of all symmetries among nfj, we can write the 

non-trivial fusion matrices as: 
[0 1 0\ 

1 n 1 n 2 n 2 

n 2 n 3 714 

\0 n 2 n 4 n 3 J 



Y 
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N 



x 



N 



x * 



/0 



V 

(° 



1 





n 2 
ri,i 
n 3 


n 2 
n 4 



1 

n 3 
n s 
m 



n 4 

m 

n 6 



n 4 

n 7 J 
1\ 

n 7 



5 



where x, ?/ are real, and z is not real. 



The modular S matrix is of the form: 
/ 1 d\ d 2 d 2 \ 

di x y y 

d 2 y Z Z 

\d 2 y z z j 

We will work on unitary modular symbols, so d\ > 1, d 2 > 1. The argument for 
general case should have only minor changes. 

The identity N X N Y = N Y N X leads to the identities: 



(A.5) 
(A.6) 

(A.7) 

N Y N X * 
(A.8) 

(A.9) 

(A.10) 



2 i 2 i 2 

n 2 + n 3 + n 4 , 



n\ + 2n3?i4, 



1 + nin 3 + n 2 (n 5 + n 7 ) 

nin 4 + n 2 (n 6 + n 7 ) = 

nin 4 + ra 2 (r?,5 + n 6 ) = ra^ + 2n 3 n 4 . 
Nx*Ny gives no new identities. But NxNx* = Nx*Nx gives us: 
n 2 n 4 + n 4 n 6 = n 2 n 3 + n 4 n 5 , 



n 5 = n 7 , 



1 + n\ + n?. 



?tJ which implies 



Case 1: n 4 = 0. 

If n 4 = 0, then n 2 n 3 = 0. First if n 2 = 0, then 1 + riin 3 
n 3 — 1, rix — 0. It follows that rii = n 2 = = 0, n 3 = 1. This leads to n\ = 2 + rij 
which has no solutions. Secondly if n 3 = 0, then n\ = 1 + which implies 
n 6 = l,n 7 = 0. Hence n 3 = n 4 = n 5 = n 7 = 0,uq = 1. This leads to n 2 = 1, 
and n\ is arbitrary. To rule out this case, notice that the labels 1,X, X* have 
exactly the same fusion rules as the rank=3 non-self dual theory. Therefore, it is a 
pre-modular category with the same fusion rules, which is necessarily modular by 
[Brj : Suppose otherwise, then (d 2 ,z,z) would be a d 2 times (l,d 2 ,d 2 ) as vectors, 
contradicting z is not real. It follows d 2 = 1, z = u for some u 3 = 1. Comparing 
the squared lengths of row 1 and row 3 of the S matrix, we see that y 2 = d\. 
Also note that d\ = 3 + nidi. Equality of the squared lengths of row 1 and row 2 
implies x 2 + 2d\ = 3. Since x is real, this does not hold if di > 0. 

Case 2: n 4 ^ 0. 
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If n 2 = 0, then n\ = 2n 3 , 1 + nin 3 = n 2 + n\. Hence 1 + n 2 = n\ which implies 
n 4 = l,n 3 = 0. So we have n\ = n 2 = n 3 = 0,n 4 = l,n 5 = n 6 = n 7 . The 
labels 1,Y form a subcategory same as the Z 2 theory, hence = l,x 2 = 1. If 
x = —1, then y = 0, and did 2 = which is a contradiction. If x — 1, then ?/ 2 = d\. 
So d\ = d\ = 1. Using g?2 = 1 + n 3^i + 2n 5 <i2 below, we see that n 3 , n 5 , hence 
^6 = ^7 = 0. So we have n\ = n 2 = n 3 = n 5 = tiq = rij = 0, n 4 = 1, which is the 
Z 4 fusion rule. 

Suppose n 4 7^ 0, n 2 7^ 0. 

The fusion rules in Table H] gives us the following identities: 
(A.ll) d\ = 1 + mdi + 2n 2 d 2 , 

(A. 12) did 2 = n 2 di + (n 3 + n 4 )d 2 , 

(A.13) d\ = n±di + (n 5 + n 6 )d 2 , 

(A. 14) d\ = 1 + n 3 di + 2n 5 d 2 . 

Combining equations, we have 
(A.15) (n 4 - n 3 )dx + (n 6 - n 5 )d 2 = 1. 

If n 4 = n 3 , then n 6 = 1 + n 2 which implies rag = 1,77-7 = 0. Hence 77.5 = 0. By 
equation ( 1A.15I) . n 5 = n G which is a contradiction. 

If n 5 = n 6 , then n\ = 1 + n| which implies n 3 = 0, n 4 = 1. Solving all equations, 
we get rii = n 2 = n 3 = 0, n 4 = 1, n 5 = n 6 = 72,7, which is the Z 4 fusion rule. 

So we may assume from now on n 2 7^ 0, n 4 7^ 0, n 4 7^ n 3 , ^5 7^ n.6- By equation 
( 1A.8P . we have 

(A.16) n 4 (n 5 - n 6 ) = n 2 (n 4 - n 3 ). 

Hence we have 

(A. 17) d2 = -d! 



n 2 n 2 (n 4 -n 3 ) 
Plugging into (1A.15j) and simplifying, we have 



(A. 18) d\ = (n l + 2nAd 1 — — . 

n 4 - n 3 

The orthogonality of S gives us: 

(A.19) x 2 + 2y 2 = 1 + 2d 2 , 

(A.20) y 2 + 2\z\ 2 = I + dj + d 2 2 , 

(A.21) (1 + x)d x + 2yd 2 = 0, 

(A.22) ydi + (1 + 2 + z)d 2 = 0, 
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d±d 2 + (x + z + z)y = 0, 



(A.24) d\ + y 2 + z 2 + z 2 = 0. 



Note that y cannot be 0. Suppose otherwise, then x = — 1, so d 2 — 0, a 
contradiction. 

The eigenvalues of N Y are d\, J^, j^-. Their sum di + Jj- — &±pl^ = rf 1 — 
^7 + ( Jj- — is an integer. The eigenvalues of N x are d 2 , J^, ^. Their sum 
c? 2 + ^ + ^ is an integer. 

If Jj- — ^ = 0, then d\ = d\. By equation (IA.12D . +d x = n 2 + n 3 + n 4 , then di, d 2 

are integers. But the sum of the eigenvalues of Ny di — Jj- is also an integer, so 
g?i = ±1. It follows that ±1 = 11-2+113+714, but n 2 , 114 are both 7^ 0, a contradiction. 
If j- — 7^ 0, then x and subsequently all y, z + z, \z\ 2 are in Q(di, d 2 ). So all 

x, y, z + z, \z\ 2 , z 2 + z 2 are in Q(di, d 2 ). By equation IA.15[ Q(di, d 2 ) is a degree< 2 
Galois extension of Q. Therefore, the Galois group of the characteristic polynomial 
Pi(t) of Ny is either trivial or Z 2 . If it is trivial, then all eigenvalues di, Jj-, 
and (i 2 are integers. So d%,d 2 ,x,y are all integers. 

From the unitary assumption, c?i, d 2 > 1. Since ^ are integers, \x\ > d\, \y\ > 
d 2 . Equation (1A.19j) implies x = ±l,y = +d 2 . Since Jj- is an integer, |x| = d\ = l. 
Then 2yc?2 = —2 implies rf 2 = 1, contradicting j- — ^ 7^ 0. 

Therefore the Galois group of is Z 2 . Since has a pair of repeated 

roots, then pi(t) is (t — m) 2 qi(t) for some irreducible quadratics qi(t) and integer 
m or (gx(t)) 2 . Assume qi{t) = t 2 + bt + c, where 6, c are integers. Note that d\ 
has to be an irrational root of pi(t). If Pi(t) has integral roots m, then = in, 
so t/ 2 > 4 x = dij- = c implies |x| > 1. By equation (1A.19[) . y 2 > d\ implies 
x 2 < 1, hence \x\ = l,y 2 = d\. It follows from equation (1A.21j) that d\ = d\. 



By equation (1A.13j) . (n 4 — 1)^! + (n 5 + n 6 )d 2 = 0. Since n 4 > 1, it follows that 



n 4 = 1, n 5 = ii 6 = 0, contradicting n 5 7^ n 6 . 

Hence pi(t) = qi(t) 2 , and di = j-, i.e. x = d\ > 1, and y 2 < d\. So the roots 

of Pi(t) are d\,d\, j^, J^. Then c?i + ^ and ^ are both integers. By equations 

CQ2TAT231 + x = -(z + z) = ^ + 1 is an integer. On the other hand, 

+ x = xfx 2 - + 1), so x = d 2 would be a rational number s if -r 2 - + 1 7^ 0. 

Then di = y^, which is also ~ fe=l=v ^ b2 ~ 4 -, but not a rational number, hence 6 = 0, 
a contradiction. If ^j 3 - + 1 = 0, then y = — ^. Substituting this and x = d\ into 



equation ( 1A.21I) . we get d 2 = 2^§ - 1. By equation ( 1A.15I) . |gQ, hence d 2 would 



be a rational number s again, a contradiction. 

Putting everything together, we have the only desired modular S matrix. 



ON CLASSIFICATION OF MODULAR CATEGORIES 



51 



References 



[Ba] P. Bantay, The Frobenius-Schur indicator in conformal field theory. Phys. Lett. B 394 

(1997) , no. 1-2, 87-88. 

[BM] D. Belov; G. Moore, Classification of spin abelian Chern-Simons theories, |hep-t h/0505235 

[BK] B. Bakalov; A. Kirillov, Jr., Lectures on Tensor Categories and Modular Functors, Univer- 
sity Lecture Series, vol. 21, Amcr. Math. Soc., 2001. 

[Br] A. Bruguiercs, Categories premodulaires, modularisations et invariants des varietes de 
dimension 3. (French) Math. Ann. 316 (2000), no. 2, 215-236. 

[CG] A. Coste; T. Gannon, Remarks on Galois symmetry in rational conformal field theories. 
Phys. Lett. B 323 (1994), no. 3-4, 316-321. 

[CP] M. Caselle; G. Ponzano, Analyticity, modular invariance and the classification of three 
operator fusion algebras. Phys. Lett. B 242 (1990), no. 1, 52-58. 

[CZ] M. Cuntz, Integral modular data and congruences, math/0611233. 

[dBG] J. de Boer; J. Goeree, Markov traces and Hi factors in conformal field theory. Comm. 

Math. Phys. 139 (1991), no. 2, 267-304. 
[DFNSS] S. Das Sarma; M. Freedman; C. Nayak; S. H. Simon; A. Stern, Non-Abelian Anyons 

and Topological Quantum Computation, arXiv:0707.1889. 
[DW] R. Dijkgraaf; E. Wittcn, Topological gauge theories and group cohomology. Comm. Math. 

Phys. 129 (1990), no. 2, 393-429. 
[ENO] P. Etingof; D. Nikshych; V. Ostrik, On fusion categories. Ann. of Math. (2) 162 (2005), 

no. 2, 581-642. 

[EK] D. Evans; Y. Kawahigashi, Quantum symmetries on operator algebras. Oxford Mathemat- 
ical Monographs. Oxford Science Publications. The Clarendon Press, Oxford University 
Press, New York, 1998. 

[F] M. H. Freedman, P/NP, and the quantum field computer. Proc. Natl. Acad. Sci. USA 95 

(1998) , no. 1, 98-101. 

[FKLW] M. Freedman; A. Kitaev; M. Larsen; Z. Wang, Topological quantum computation. Bull. 

Amer. Math. Soc. (N.S.) 40 (2003), no. 1, 31-38. 
[FKW] M. H. Freedman; A. Kitaev; Z. Wang, Simulation of topological field theories by quantum 

computers. Comm. Math. Phys. 227 (2002), no. 3, 587-603. 
[FLW1] M. H. Freedman; M. J. Larsen; Z. Wang, A modular functor which is universal for 

quantum computation. Comm. Math. Phys. 227 (2002), no. 3, 605-622. 
[FLW2] M. H. Freedman; M. J. Larsen; Z. Wang, The two- eigenvalue problem and density of 

Jones representation of braid groups. Comm. Math. Phys. 228 (2002), 177-199. 
[FW] M. H. Freedman; Z. Wang, Large quantum Fourier transforms are never exactly realized 

by braiding conformal blocks. Phys. Rev. A (3) 75 (2007), no. 3, 032322, 5 pp. 
[FQ] D. Freed; F. Quinn, Chern-Simons theory with finite gauge group. Comm. Math. Phys. 156 

(1993), no. 3, 435-472. 

[FK] J. Frohlich; T. Kerler, Quantum groups, quantum categories and quantum field theory. 

Lecture Notes in Mathematics, 1542. Springer- Verlag, Berlin, 1993. 
[FNTW] E. Fradkin; C. Nayak; A. Tsvelik; F. Wilczek, A Chern-Simons Effective 

Field Theory for the Pf affian Quantum Hall State, Nucl.Phys. B516 (1998) 704-718, 

|arXiv:cond-mat/9711087| 
[FTL] A. Feiguin; S. Trebst; A. W. W. Ludwig; M. Troyer; A. Kitaev; Z. Wang, M. Freedman, 

Interacting anyons in topological quantum liquids: The golden chain, Phys. Rev. Lett. 98, 

160409 (2007). 

[G] T. Gannon, Modular data: the algebraic combinatorics of conformal field theory. J. Alge- 
braic Combin. 22 (2005), no. 2, 211-250. 



52 



ERIC ROWELL, RICHARD STONG, AND ZHENGHAN WANG 



[GK] D. Gepner; A. Kapustin, On the classification of fusion rings. Phys. Lett. B 349 (1995), 
no. 1-2, 71-75. 

[HH] T. Hagge; S. Hong, Some non-braided fusion categories of rank 3, arXiv: 0704.0208 
[HRW] S.-M. Hong; E. C. Rowell; Z. Wang, On exotic modular tensor categories, Commun. 

Contemp. Math. Vol. 10 (2008), Suppl.l, 1049-1074, arxiv: 0710.5761. 
[Kil] A. Kitaev, Fault-tolerant quantum computation by anyons. Ann. Physics 303 (2003), no. 

1, 2-30. 

[Ki2] A. Kitaev, Anyons in an exactly solved model and beyond. Ann. Physics 321 (2006), no. 1, 
2-111. 

[KL] L. Kauffman; S. Lins, Temperley-Lieb recoupling theory and invariants of 3-manifolds. 
Annals of Mathematics Studies, 134. Princeton University Press, Princeton, NJ, 1994. 
x+296 pp. 

[KW] D. Kazhdan; H. Wenzl, Reconstructing monoidal categories. I. M. Gelfand Seminar, 1 11 

136, Adv. Soviet Math., 16, Part 2, Amer. Math. Soc, Providence, RI, 1993. 
[LRW] M. J. Larsen; E. C. Rowell; Z. Wang, The N -eigenvalue problem and two applications, 

Int. Math. Res. Not. (2005), no. 64, 3987-4018. 
[LWa] M. Larsen; Z. Wang, Density of the SO (3) TQFT representation of mapping class groups. 

Comm. Math. Phys. 260 (2005), no. 3, 641-658. 
[LWe] M. Levin; X.-G. Wen, String-net condensation: A physical mechanism for topological 

phases, Phys. Rev. B71, 045110 (2005). |cond-mat/040461~7| 
[Ma] M. Manoliu, Abelian Chern-Simons theory. I. A topological quantum field theory. J. Math. 

Phys. 39 (1998), no. 1, 170-206. 
[Ml] M. Miiger, From subfactor to categories and topology, II J. Pure Appl. Algebra 180 (2003), 

no. 1-2, 159-219. 

[M2] M. Miiger, On the structure of modular categories, Proc. London Math. Soc. (3) 87 (2003), 
no. 2, 291-308. 

[MR] G. Moore; N. Read, Nonabelions in the fractional quantum hall effect, Nuclear Physics B, 

Volume 360, Issue 2-3, p. 362-396. 
[MSI] G. Moore; N. Seiberg, Lectures on RCFT. Superstrings '89 (Trieste, 1989), 1-129, World 

Sci. Publ., River Edge, NJ, 1990. 
[MS2] G. Moore; N. Seiberg, Classical and quantum conformal field theory. Comm. Math. Phys. 

123 (1989), no. 2, 177-254. 
[NR] D. Naidu; E. C. Rowell, A finiteness property for braided fusion categories, arXiv:0903.4157, 
[01] V. Ostrik, Fusion categories of rank 2. Math. Res. Lett. 10 (2003), no. 2-3, 177-183. 
[02] V. Ostrik, Pre-modu lar categories of rank 3, Mosc. Math. J. 8 (2008), no.l, 111—118, 

|arXiv:math/0503564| 

[P] J. Preskill, Chapter 9 at |http: / / www.theory.caltech.edu/ ^pr eskill / ph229 / 

[Rol] E. C. Rowell From quantum groups to unitary modular tensor categories in Contemp. 

Math. 413 (2006), 215-230. 
[Tu] V. Turaev, Quantum Invariants of Knots and 3-Manifolds, De Gruyter Studies in Mathe- 
matics, Walter de Gruyter (July 1994). 
[TV] V. Turaev; O. Viro, State sum invariants of 3-manifolds and quantum 6 j -symbols. Topology 

31 (1992), no. 4, 865-902. 
[TW] I. Tuba; H. Wenzl, On braided tensor categories of type BCD J. reine angew. Math. 581 

(2005), 31-69. 

[Wa] Z. Wang, Topologization of electron liquids with Chern-Simons theory and quantum com- 
putation. Differential geometry and physics, 106-120, Nankai Tracts Math., 10, World Sci. 
Publ., Hackensack, NJ, 2006. 



ON CLASSIFICATION OF MODULAR CATEGORIES 



53 



[Wal] K. Walker, On Witten's 3-manifold Invariants, 1991 notes at http://canyon23.net/math/ 
[Wen] X.-G. Wen, Topological Orders and Edge Excitations in FQH States, Advances in Physics 

44, 405 (1995). |cond-mat/9506066[ 
[WW1] X.-G. Wen; Z. Wang, A classification of symmetric polynomials of infinite variables- a 

construction of Abelian and non- Abelian quantum Hall states, Phys. Rev. B, 77, 235108 

(2008); cond-mat/0801.329 
[WW2] X.-G. Wen; Z. Wang, Topological properties of Abelian and non- Abelian quantum Hall 

states from the pattern of zeros, Phys. Rev. B, 78, 155109 (2008). larXiv:0803.1016l 
[Wenz] H. Wenzl, C* tensor categories from quantum groups. J. Amer. Math. Soc. 11 (1998), 

no. 2, 261-282. 

[Wil] F. Wilczek, Fractional Statistics and Anyon Superconductivity, World Scientific Pub. Co. 
Inc. (December 1990). 

[Witt] E. Witten, The search for higher symmetry in string theory. Physics and mathematics of 
strings. Philos. Trans. Roy. Soc. London Ser. A 329 (1989), no. 1605, 349-357. 
E-mail address: rowell@math.tamu.edu 

Department of Mathematics, Texas A&M University, College Station, TX 77843, 
U.S.A. 

E-mail address: stong@ccrwest.org 

Center for Communications Research, 4320 Westerra Court, San Diego, CA 
92121-1969, U.S.A. 

E-mail address: zhenghwa@microsoft.com 



Microsoft Station Q, CNSI Bldg Rm 2237, University of California, Santa 
Barbara, CA 93106-6105, U.S.A. 



